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STATISTICAL CORRELATION ANALYSIS FOR COMPARING VIBRATION DATA 

FROM 

TEST AND ANALYSIS 


INTRODUCTION 

Up until the last 30 years, structural analysis consisted largely of abbreviated and approxi- 
mate calculations. Typically, the calculations represented empirical knowledge more than an un- 
derstanding of structural behavior. Limited analytical capabilities also lead to limitations in test- 
ing. Without the knowledge of how a structure should respond to loads, the testing of a new de- 
sign often consisted simply of putting it to its intended use. If it worked sufficiently well for a 
sufficiently long time, it was deemed successful and became a model for future structures of a 
similar nature. 

This situation changed drastically as the result of three major developments. With the advent 
of the digital computer in the 1940s came the possibility of doing the large number of calculations 
required to determine accurately and in detail the stresses and deformations that a structure would 
experience when serving its intended functions. Next, the development of the finite element me- 
thod in the 1950s provided algorithms, executable on digital computers, which could indeed calcu- 
late stresses and deformations in complex structures. Finally, software systems were developed to 
combine structural theory, finite element methods, and computer hardware into useful tools for 
structural analysis. 

The NASTRAN software system (NASA STRuctural ANalysis), development of which was 
begun by NASA in 1965, is an example of such a system which has received wide usage. There 
were a number of reasons for its acceptance. First, it integrated into a unified set of software a 
sufficient number of algorithms and data types to carry out the majority of structural analyses. 
Second, NASTRAN included software to overcome many of the deficiencies of early computers, 
such as limited memory size, which would have otherwise limited the complexity of finite element 
models. This meant that the limitation on the complexity of an analyst’s model shifted to the 
amount of computer time he could obtain. Finally, extensive effort went into providing good user 
interfaces, complete documentation, and wide dissemination. To enable wide use, NASTRAN is 
written mainly in a computer independent subset of FORTRAN, which permits it to run on 5 dif- 
ferent makes of computers. The total capabilities provided by NASTRAN have contributed sig- 
nificantly to the present situation, where finite element analysis using digital computers is not only 
the dominant means of accomplishing structural analysis, but it is also widely used for analysis in 
heat transfer, aerodynamics, acoustics, control systems, potential theory, and actually any prob- 
lem which can be solved by sets of linear equations. 

The overall effect of this new tool has been to improve the quality of structural designs. With 
Finite element techniques, the designer can be much more certain of reaching a desired tradeoff of 
structural efficiency, cost, and reliability than previously. However, the volumes of data produced 
by complex finite element analyses, now routinely performed, have created problems of inter- 
pretation, because of the human time required to assimilate large amounts of data. For the same 
reason, test data is often not completely compared to predicted analytical results. Often, if a 
structure goes through static and dynamic testing with no apparent problems, such as breakage or 
pronounced deflection or vibration, it is deemed adequate. When test data obviously disagrees 
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with results predicted by finite element analysis, a too frequent practice is to assume that the ana- 
lytical data is wrong and to adjust its parameters until satisfactory agreement with test data is 
achieved. 

There are a number of ways to increase the value of combined analytical and test data on 
structural behaviour. The work described here was undertaken to derive some of these improve- 
ments. Fundamentally, the two sets of data provide a means of cross checking, thereby reducing 
the probability of undetected errors in either set. A principal problem in doing cross checking is 
that there is not a commonly accepted measure of comparison between analytical and test data. 
Often, this means that a few hundred data points from both sources, regarded as critical, are plot- 
ted in some way, and if the data plots appear reasonably alike to the eye, or can be made to look 
alike by reasonable alterations to the analytical model, then agreement is decreed. 

With the techniques described below, it is possible to use NASTRAN to place analytical re- 
sults in a format in which they can be statistically compared as a group to equivalent groups of test 
data gathered by automated test systems. Statistical routines can then provide criteria describing 
overall agreement and identify areas of significant disagreement. This means that far more data 
can be evaluated with far less human labor. This technique can help identify whether disagree- 
ment is due to faulty test procedures or faulty modeling. However, another promising long range 
benefit of the statistical correlation of analytical and test data is that structural design can be fur- 
ther improved. There is presently no widely accepted and effective way of validating analytical re- 
sults on real problems. Statistical comparisons with equivalent test data appear to provide an ac- 
ceptable means for this validation. If this technique becomes widely used and improved, then it 
will gradually lead to improved analytical techniques, thus leading in turn to improved designs and 
eventually to improved structures. 


OBJECTIVE 


When designs need to qualify under narrow allowances, as is normally the case with aerospace 
structures, it is often necessary to compare structural vibration properties obtained from testing 
with those obtained from analysis. 

In order for a comparison to be meaningful, it is necessary to define the basis on which the 
data will be compared. In general, they will be compared on the basis of like modes. A test mode 
and an analytical mode are alike if the amplitude of the standing wave pattern is “exactly” the 
same for both at all sample points. The bothersome qualifier in this definition of like modes is 
the word ‘exactly’. Even if one were to assign a reasonable tolerance for defining the word ‘exact’, 
it would be rare to find a test mode that measured exactly the same as an analytical mode. A 
more realistic definition of modes being “alike” is to require that the pair correlate closely. This 
definition allows some engineering judgement for prescribing the tolerance of “closeness.” 

There are no exact solutions to the vibration characteristics of complex three dimensional 
structures. Highly capable analytical techniques are available, but they are liable to human error. 
Highly capable experimental equipment is also available, but human error can enter the experimen- 
tal side from two sources: in fabricating the test article and in conducting the test. Once a like 
pair is found, it can be expected that the frequency at which this mode appears is different for 
analysis than it is for test. Therefore, the purpose of this enterprise is not to judge the correctness 
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of either set of data to represent a given design, but to point up differences and equip the respon- 
sible structural engineers with enough measures for them to be able to assess the areas of disagree- 
ment. 

To illustrate how statistical comparison can help the structural engineer, consider the follow- 
ing example. After testing a collection of modes for likeness in a given frequency range, say 8 
from test and 1 1 from analysis, it might be that 7 matches are found. Also, consider that one test 
mode and four analytical modes are not matched. In those that match, it is possible that the fre- 
quencies of the analytical modes are all higher than their test counterparts. It may also be that the 
one test mode without a match measures higher in frequency than all other test modes, while two 
of the analytical modes without matches measure lower in frequency than all other analytical 
modes, and the other two are in the middle range. Figure 1 depicts this example. 


Analysis 



Figure 1 

One would be inclined to look above the upper frequency bound of the analytical modes to Find a 
match for the spinster test frequency, and to look below the lower frequency bound of the test 
modes to find matches for the mateless analytical modes. As for the two unpaired analytical modes 
in the middle frequency range, there is a possibility that analysis uncovered two modes that test 
failed to find. It might also be that one grid point stands out in five of the matching pairs of modes 
as the place where the modes show the greatest difference in their shapes. One might be suspicious 
of the analysis as having too low a value of mass at that point. 

The happy ending to this story might be that the analytical mass defect was found, and the 
wanting test instrumentation was supplied so that the overall shift in frequency was erased, two 
middle range test modes were revealed, and 1 0 out of 1 1 matches were found in the frequency range 
on a second trial. 

The limited amount of experimental data and the limited access in which to install sensors at 
interior points are factors that constrain the points at which comparisons can be made. There is a 
disproportionately large number of points at which analytical data is available for comparisons. 
The objective will be to bring about useful comparisons within the constraints of available data. 
The intent is to be able to identify the mode shapes from both of the sources; to say how closely 
any two mode shapes are related; to say at what frequencies the test mode and the analytical mode 


3 



of a closely related pair appear; to give the locations at which two modes diverge beyond a thresh- 
old; and to isolate the greatest divergence. 

A correlary objective is to establish a communication basis between the records from analysis 
and the records from test. Once such a channel is established, then standardized data can be ex- 
changed by means of a program. For instance, the arrays of generalized mass and generalized stiff- 
ness can be passed from analysis to test while the array of modal damping can be passed from test 
to analysis. 


METHOD 


Structural properties determined by either analysis or test can properly be regarded as non- 
deterministic and therefore random. For instance, in analysis the factors liable to random varia- 
tion are: the interpretation of drawings, the modeling of elastic relations, the representation of 

joints, the modeling of mass properties, entering of data into the computer, transfer of data into, 
out of, and within a computer, computer behavior, and debugging methods. There will be as many 
analytical results as there are analysts. Similarly, variability in testing can arise from the quality 
of materials in fabrication, the machining of parts, the assembly of components, the mounting of 
sensors, the calibration of recorders, the reading of records, and the processing of data. The ex- 
perimental results will be as individual as the number of teams conducting tests. Since it is the 
intention to compare two random records, a well proven stochastic tool to apply is the correla- 
tion function— or particularly-the correlation coefficient. Thus, the criterion for deciding wheth- 
er two modes are the same can be mediated by how close to unity the correlation coefficient for 
them is. Statistical methods will be employed to compare test mode shapes with analytical mode 
shapes. Analytical and test measurements of mode shapes will be taken at the same set of struc- 
tural locations, so that behavior will be compared at like points. At present, the comparisons in 
this study will be limited to undamped real modes, as opposed to damped complex modes. A test 
mode shape and an analytical mode shape will each be considered to be a statistical set of data. 
The number of test modes to be compared need not be the same as the number of analytical 
modes, so that all available modes from both sources within a specified frequency range can be 
considered. The primary tool for making comparisons will be the correlation coefficient. 

According to James & James ‘Mathematical Dictionary’, the definitions of the statistical 
quantities that will be used in this analysis are as follows: 

n 

Variance = 2 (xj — x')^ p(xj) = (1) 

i=l 

where n is the number of sample points in the statistical file, 

1 n 

x' is the mean = — X x; (2) 

n i=l 

and p(Xj) is the probability of finding the sample Xj in any statistical record. Other statistical 
quantities will be introduced after discussing those introduced thus far. 
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The following interpretations will be given to these quantities for purposes of this analysis. 
Data from each mode is a distinct statistical record. The quantity Xj is the measurement of the de- 
parture from the equilibrium position at a given point “i” of a single vibration mode from an eigen- 
value analysis. If xj were measured in time from an undamped oscillation, it is known from the 
definition of a mode shape that there will be an equal number of positive excursions as negative 
excursions away from the equilibrium position of a mode. The mean of a time varying sample 
would be the equilibrium position, so x' would be zero. The Xj’s are not taken as time varying 
samples, but as the normalized envelopes of a finite number of extreme model positions for 
which the average is definitely not zero. The probability p(xj) is a quantity which varies from 
mode to mode. The function p(xj) for a mode is not known a priori, so it will be estimated by 
using the scheme of empirical probability density as described in Y. W. Lee’s “Statistical Theory of 
Communications”, page 125. The range of Xj depends on how a mode is normalized. For example, 
if a mode is normalized to its maximum component, the range is — 1 < xj < + 1. For normaliza- 
tion according to a particular component or to mass, the ranges can be greater or less than unity. 
The range is divided into equal increments, Ax. Modal data is scanned for how many of the modal 
displacements fall within each increment; i.e, the number Nj in each increment Axj. The frequency 
ratio is constructed for each increment, by dividing the count in each increment Nj by the total 

number of modal data points n: Nj/n. According to Jacob Bernoulli’s theorem on probability, 
N,/n 

the quantity 1 stays within an e neighborhood of p(xj) as n gets large and Ax grows small. It 

Ax 

will be assumed that the number of samples of a mode shape n is sufficiently large to justify the 
Nj/n 

use of the quantity as a reasonable approximation to the probability density function p(x) 

Ax 

for use in this algorithm. The probability that a model sample point “h” will fall within the 
range of 

x + Ax 

x<h<x + Axis J*P(x)dx 

X 

Because of the discretized nature of the available data the integral will be approximated in incre- 
mental form by 

P(Xj < h< xj + Ax) as (xj + Ax - xj) ** Nj/n (3) 

Ax 

On occasion, a fault in a sensor or a recorder will arise which voids data for some point. This ren- 
ders the probability null for finding data at that point. Instead of admitting a zero probability, 
this analysis will reduce the sample space “n” to the corresponding number of data points avail- 
able in both analysis and test. 

With the above provisions for this analysis, the statistical quantities reduce to: 

Variance 

n 

° 2 = 2 ( x i - x 1 ) 2 Nj/n (4) 

i=l 
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Standard Deviation 


a = +\^2 = + [ 2 (Xi-x') 2 Ni/n]' /2 

i=l 


(5) 


R.M.S (Root Mean Square) 


P = U/n£ X} 2 Nj/n] 1//2 . If x' = 0 then p = a. 
Correlation Coefficient between two statisitical files xj(a) and xj(e) is: 


*ab 


2 (Xj(a) - x' a ] P a (x i )][x i (e) - x' e ] P e ( x i> • Pae( x a> x e)i 
= JSl — 




i(a) - x' a ] 2 p a (xj) 2 1 x i^ e ) - x 'el Z M x i) 
i=l 

using x' a = x' e = 0, and using the frequency ratios for approximating probability causes “r” to sim- 
plify to: 

N a (xj) N e (xj) 1 

2 j Xj( a ) Xj(e) -y N a ( Xi ) N e ( Xi ) 


l 2 


2 Xj(a) xj(e) 


i=l 


n n 


i=l 




*ab 


2 Xj(e) 2 i!£2 

i=l n i=i n 


a a °e 


n 

r ab = 2 ' 
i=l 


Xj(a) 

°a 


Xj(e) 

°e 


N a (xj) Ne(xj) 


( 6 ) 


Equation (6) can be interpreted as the total effect over the whole mode of multiplying the 
modal amplitude in the “a” set of data when scaled by its standard deviation, by the modal am- 
plitude in the “e” set of data when scaled by the “e” standard deviation, for every pair of corre- 
sponding points, and summed for all n points. This method gives complete freedom to both the 
analyst and to the experimenter for organizing data in the manner of his choice (so long as it is 
uniform) and of normalizing to any convenient scheme independent of each other. When each set 
of data is scaled to its own standard deviation, both are placed on equal footings regardless of how 
the modal data was originally obtained. 

All modes in both sets will be compared with each other, and the correlation coefficient will 
be computed for every combination. 
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Relative Deviation 


It is instructive to determine where and by how much test and analytical mode shapes differ 
from one another. The approach here is to scale the modes so as to be highly sensitive to differ- 
ences. Scaling by the RMS value for each mode will put both modes on such a footing that their 
differences will show a greater spread than they would if scaled by the variance. Using subscript 
“a” to indicate quantities belonging to analytical data and “e” to indicate quantities belonging to 
experimental data, the expression for relative deviation in the i^ 1 sample of the j^ 1 mode is writ- 
ten as: 


(Rel. Dev.)jj 


xi(a)j 


P 


a 


Xi(e)j 


(7) 


Relative deviation in a mode is the absolute value of the difference between the analytical 
modal displacement scaled by its modal RMS, and the experimental modal displacement scaled by 
its own modal RMS. An option will be given to specify what percentage of spread is to be dis- 
played with the grid point identity. The default value of this spread is 5%. 

Three other quantities will be computed which have been found to be useful by practicing 
test engineers. These quantities are defined in an internal GSFC memorandum written by D. J. 
Hershfeld dated January 8, 1971, entitled “Fitting Mode Shapes to Experimental Data.” Accord- 
ing to equation (2) of the memo (when transcribed into notation adapted in this report) the nor- 
malizing factor of the j^ 1 mode is: 


Normalizing Factor 
n 

E x i(a)jXi(e)j 

q (8) 

J n 

X Xj(a)j 2 
i=l 

The standard error, according to equation (5) of the memo in transcribed notation, is: 

Standard Error 

n n 

S = | - [ 2 xj(e) 2 N e (Xj) - 2 x^a) 2 N a ( x >) ] | % = I a e 2 - q 2 | % (9) 

" i=l i=l 

Differences in mode shapes are computed by a common scaling of their absolute values by the 
standard error: 


Scaled Difference = | [xj(a) — Xj(e)]/S | 


( 10 ) 
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The tabulations of input modal arrays will be processed by this program. Listings will be made 
of modal masses, modal stiffnesses, and modal damping arrays versus frequency. 


IMPLEMENTATION 


In order to put this statistical correlation method of comparing test and analysis into op- 
eration, attention must be given to three distinct operations: a. the generation of the analyt- 

ical vibration data and the transfer of its results, b. the sensing of the test vibration data and 
the transfer of its results, and c. the computation of the comparison between test and analysis 
and the acceptance of data from the two sources. 

Analytical vibration data are expected to be implemented with NASTRAN. Output of the 
results is needed in non-standard form, so an ALTER packet was written to tailor the data and 
to translate the data to BCD for universal readability. These analytical data need further pro- 
cessing before computation can start and are passed to intermediate processors. 

Two assemblies of test data are expected. They are: the list of points where test data 

were taken and the tabulation of mode shapes and frequencies in simple BCD format. They 
need to be transferred to files in the primary computer for computation. 

The computation according to the theory outlined above is accomplished in a program 
called STATCORR. It expects input data in particular formats and allows the user to choose 
from a variety of output options. 

Consideration has been given to analytical and test data coming from various sources. The 
STATCORR program will operate on data independent of its source, so long as it is entered ac- 
cording to a specified format. During this development, the analysis was done with NASTRAN 
and the tests were performed by the Environmental Test & Integration Branch at GSFC. Refer- 
ences to specific analytical operations will pertain to capabilities in NASTRAN. 

Provision was made for data coming from several runs, structures with symmetry, different 
normalizations, and non-uniform samples. 

Liaison 


At the outset, it is mandatory that comparisons be based on data gathered from identical 
locations. Liaison between test and analysis personnel prior to either the generation of the NAS- 
TRAN model or the instrumentation of the test article must establish the points and the compo- 
nent directions at which data is to be sampled. This set of locations will be referred to as the 
TESET vector. TESET will govern the location and orientation of test sensors and will govern 
the entries in the analytical partitioning vector. It is used in tabulating mode shapes from test 
and analysis, and it is the primary reference within STATCORR for sorting and reporting differ- 
ences. 

Equally important as matching data points is matching boundary conditions. The most dif- 
ficult experimental boundary to enforce is clamped. The easiest to enact is free-free. Liaison 
between the test engineer and the analyst will help to uncover any disparity between the test 
set-up and the analytical boundary constraints. 
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Provision has been made for using symmetry. If the structure has symmetry that can be 
taken advantage of, liaison should define whether both test and analysis will collect data on only 
the basic segment. If so, agreement is needed on the naming of each plane of symmetry. Then, 
each test mode and each analytical mode must be tagged with SYM # PLAN parameters for the 
type of symmetry across each plane of symmetry. If only analytical data is based on symmetry, 
while test data is gathered over the whole structure, liaison is needed for three additional items. 
Octants of symmetry need to be identified. Each test point must be checked for the existence of 
an image point of reflection in the basic segment of analysis. Multiply-reflected image points must 
be provided with clones. A list of points with their individual component numbers belonging to 
each octant is organized by analysis, but should be double checked by the test engineer before 
proceeding. 

Analysis 


In NASTRAN, there are 5 options for eigenvalue analysis: Determinant, Inverse Power, 

Givens, Feer, and Hessenberg. All methods are manageable for use in conjunction with STAT- 
CORR. There is a problem, however, in transferring NASTRAN eigenvalue results into STAT- 
CORR. A DMAP ALTER has been written that manages this transfer. Modal frequencies, mass, 
and stiffness are delivered as part of the LAMA table in DTI format, using TABPCH. Parti- 
tioned mode shapes are delivered as a single matrix in DMI format using OUTPUT3. An explana- 
tion of the ALTER packet will follow subsequently. There are some precautions that should be 
observed with some of the 5 methods. With the Feer method one should check-point the eigen- 
value run and exit after the READ module. Results should be examined for the acceptability of 
the last half of the roots extracted. If negative roots or vectors that fail to meet the orthogonal- 
ity criterion are found, the PHIA matrix and LAMA table should be partitioned down to only 
the valid roots immediately after restarting; then one can proceed with the prepared DMAP 
ALTER route. Sometimes with the Determinant and the Inverse Power methods, several runs may 
be required before all the roots in a given frequency range are found. Data from any number of 
separate runs may be combined, but some roots may be repeated in different runs. A post-proces- 
sor program (MERGE) has been provided to combine eigenvector data from the several runs, de- 
tect and eliminate duplicates, and sequence them in ascending order of frequencies. Modes may 
be normalized according to MASS, MAX, or POINT, but once a decision is made for a method of 
normalizing, it must be used exclusively for all analytical modes henceforth. 

Partitioning is done using the DMAP module PARTN which requires a user supplied parti- 
tioning vector. The vector has been arbitrarily named TESET (standing for test set), but this 
name is imbedded in the DMAP statements of the ALTER packet, so it becomes necessary to use 
this name. DMI bulk data cards define the vector as a rectangular matrix of a single column, and 
the number of rows equal to order G of the structural model. This partitioning vector will be 
designed to serve a double duty. In addition to its duty of partitioning the PHIG matrix, TESET 
will serve the STATCORR program as a guide for correlating the sequence of test data with ana- 
lytical data. The device employed to enable this secondary duty is that of constructing the DMI 
vector in ordered pairs. The first member of this pair is a NASTRAN internal sequence number 
which corresponds to a component of a test point. The second member of the pair is the Grid 
Point ID number belonging to the test point. The internal sequence number is actually the row 
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number of the equation of equilibrium for a given degree of freedom as they are organized in 
NASTRAN’s system of matrices. The matrix of eigenvectors, PHIG, is the matrix to be parti- 
tioned, so the vector used to do the partitioning has to be the same order as the subject matrix. 
The final letter of standard NASTRAN matrix names gives its size; so in this case, the partition- 
ing vector must be of order G; or briefly, the vector is G-sized. The format of the DMI input 
data requires that the first member of the pair is to be an integer, and second member of the pair 
is to be a real number. Consequently, in implementing this DMAP ALTER the user is required to 
first list the grid point identification number (ID) as a real number. Just placing a decimal point 
after the ID number makes it a real number. The internal sequencing used by NASTRAN for 
matrix operations could be different than that which the analyst originated with his grid point 
sequencing. The partitioning vector must be organized according to internal sequencing, so it 
is advisable to include DIAG 21 and 22 and CHKPNT YES in the executive control packet, and 
then insert an ALTER to EXIT somewhere after the GP4 module. The detailed correspondence 
between internal and external sequencing for the various sets will be printed when DIAG 21 and 
22 are activated. The internal sequence number of every corresponding entry on the DOF list 
will become the row number of the partitioning vector. The row numbers are required to appear 
in ascending sequence on the DMI cards. 

Symmetry 


As a special analytical case, an analyst likes to take advantage of symmetry of a particular 
structure for the economy and improved definition of modes. Cyclic symmetry is automated 
within NASTRAN, but reflective symmetry is not. Mode shapes analyzed via cyclic symmetry 
can automatically be described either for the whole structure, or for just a segment. The only 
additional processing needed for such modes, before they would be ready for STATCORR, is 
partitioning. On the other hand, modes analyzed via reflective symmetry require additional 
processing in order to prepare the modes for STATCORR. Eigenvalue analyses employing re- 
flective symmetry may cause the analytical data to be delivered in the form of several runs, de- 
pending on the available combinations of symmetry and anti-symmetry. The post processor pro- 
gram for analytical data, mentioned above in the discussion of DET and INV methods, will operate, 
as explained before, by combining all modes, eliminating duplicates, and sequencing both frequency 
and GP list. In addition, the analyst will supply 3 parameters that will act like flags to the STAT- 
CORR program for indicating the symmetry status with respect to each of the bounding planes of 
symmetry. Each of the three parameters can carry the values — 1 , 0, or +1 for the symmetry status. 
“0” denotes ‘'absence of symmetry” such as for the face plane of a plate. Positive 1 denotes sym- 
metric behavior about the plane of symmetry. Negative 1 denotes anti-symmetric behavior about 
the plane of symmetry. Names given to these three parameters are SYM1PLAN, SYM2PLAN, and 
SYM3PLAN. Values for these symmetry parameters must always be supplied, if modal data is 
based upon information gathered on the basis of only a basic segment of a structure. If data from 
both test and analysis are flagged for OK having symmetry, STATCORR tests for the existence of 
symmetry for one file and anti-symmetry for the other file occurring simultaneously with respect 
to any one of the three planes. When such an unlikeness occurs, the algorithm immediately de- 
clares the correlation coefficient to have the value zero, and suspends calculations for both vari- 
ance and RMS. Likeness amongst all three parameters between the two sets simultaneously indi- 
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cates to STATCORR that only the correlation coefficient for the basic segment need be calcu- 
lated, because it is equal to that for the complete structure. RMS is also the same for the whole 
structure as it is for the basic segment. However, variance of the whole structure is equal to the 
number of octants of symmetry multiplied by the variance of the basic segment. Standard devi- 
ation of the whole structure is equal to the square root of the number of octants of symmetry 
multiplied by the standard deviation of the basic segment. 

To illustrate how the correlation coefficient has values alternating between zero and that for 
the basic segment, a pair of free parallelepipeds with 3 planes of symmetry will be used. Modes 


SYM 



involving translations in only one coordinate will provide enough detail to be realistic without 
obscuring the concepts with complications. 

Symmetry parameters SYM 1 PLAN, SYM2PLAN, and SYM3PLAN for the two solids being 
compared have the values +1, —1, +1 for solid A and +1, +1, +1 forB. A set of modes conforming 
to these parameters is diagrammed to create a mental image. 

The correlation coefficient between the two structures is 
n 

X x m 


i=l 

IX x i 2 l ,/2 l X y i 2] ^ 

i=l i=l 


(xy)i + (xy) 2 + (xy)3 + (xy>4 (xy)s 

[x] + x 2 + x| + . . . . \\}' A [y] + y\ + y\ + . . . . y \Y A 


n s n s _ n s 

where £ Xjyj = (XY) S , £ xj 2 = X s , and £ yj 2 = Y s - 

i=l i=l i=l 

Substitute the equivalent, signed basic segment portion for octants 2 and above. 

(xy)] + (xy)i + [( — x)y] ; + [( -x)y] ; + (xy)j + (xy)i + [( -x)y] ] + [( -x)y] i 

[ 8x ! ] 1/2 [ 8 Y1 ] 1/2 

0 

0. 

8a x o y 
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If, however, both A and B were symmetric with respect to all 3 planes of symmetry, the correla- 
tion coefficient can be represented in terms of basic segment data as follows: 


8(xy)j 


8(xy) j (xy) j 


r = 


8 1/2 o x 8 1/2 Oy 


8 C7 X Oy 


q.e.d. 


( 11 ) 


a x Oy 


This concept is built into the correlation algorithm of STATCORR. A table of correlation coeffic- 
ients for various symmetry permutations clearly shows the oscillation between two values. 


SYMMETRY 

PARAMETERS 


CORRELATION 

COEFFICIENT 


SYMMETRY 

PARAMETERS 


CORRELATION 

COEFFICIENT 


S 0 

0 : x 

0 

S A 0 : 

X 

0 

A 0 

O : y 

A A 0 : 

X 

A 0 
S 0 

0 : x 

0 : x 

0 

A A 0 : 
A A O : 

X 

y 

BASIC 

A 0 
A O 

O : x 
0 : y 

BASIC 

S A S : 
S S S : 

X 

y 

0 

S O 
S 0 

O : x 
O : y 

BASIC 

S S S : 
S S S : 

X 

y 

BASIC 

S A 

0 : x 

0 

ASS: 

X 

0 

S S 

0 : y 

S A A : 

y 

S S 
S S 

0 : x 
0 ; y 

BASIC 

ASA: 

ASA: 

X 

y 

BASIC 

A S 

0 : x 

0 

S A S : 

X 

0 

A A 

0 : y 

S A A : 

y 

A S 
A S 

0 : x 
0 : y 

BASIC 

AAA: 

AAA: 

X 

y 

BASIC 

Variance. The variance is given by 





n 






- 2 X|2 






i=l 






a 2 = Xj + X 2 + X 3 + X 4 + X 5 + X 6 + X 7 + X 8 
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which, in terms of equivalent basic segment for the parallepiped, becomes: 


a 2 = 8 Xj = 8 a, 2 , or a 2 who]e = 2? a 2 basic 


where p = the number of planes of reflective symmetry. 


( 12 ) 


Standard Deviation. Similarly, 


° - ( ff2)^ = (8 a , 2)^ = 2\ZI a j or a,, = 2P/2 


whole ^ ' " a basic 


Root Mean Square 


p = ojyj n , where m equals the sample points in the basic segment and n= 8 m. 
p basic ” a basic V® > but a whole _ 2\/5" ^asic 


p whole V* a b as i c /V^ a basic p basic 

p whole ~ p basic 


( 13 ) 


Normalization Factor 


C = 


Same numerator as “r” 


Using the results leading up to equation (11), one can conclude that when there is unlikeness in 
the symmetry parameters, the numerator is zero, causing C to be null When there is likeness 
amongst all of the symmetry parameters, 


C = 


8(XY)i _ (XY), 


8Xj 


Xl 


( 14 ) 


C - c 

whole basic 


Standard Error 


s= [L(o e 2 -o a 2 )] % 


13 



then 


a e 2 - 8Xj(e) and a a ^ = 8Xj(a), Make these substitutions. 


S = [— (X 1 (e)-X 1 (a)] 1/2 
8m 


^ m ^basic^ 


*/2 


^whole ^basic 

The final alternate of reflective symmetry is the case in which analytical properties are determined 
by only modeling a basic segment, while test uses the entire structure. The scheme will be to limit 
the calculations to only those quantities that are necessary. It was discovered during the develop- 
ment of this method that it was not necessary to deploy the basic segment geometry into a mode 
shape over all 8 octants. By planning ahead, it is possible to arrange that every instrumented point 
on the test article, when reflected back to the basic segment, will define an image point in the 
basic segment which is a member of the grid of points used in the analysis. It is possible that sev- 
eral external test points reflect back to a common image point in the basic segment, but in the 
operation of the processor program, it is necessary to provide a unique point in the basic segment 
for every instrumented point. A multiply-reflected point violates this requirement. This dilemma 
is solved by setting up clones of grid points in the NASTRAN bulk data. Thus, the total set of 
points used by analysis in the basic segment must embrace all of the instrumented points. Those 
will be of three kinds: direct, image, and cloned points. The direct points correspond to the 

instrumented points in the basic segment on the test article The image points are those having no 
corresponding points in the basic segment, and that, when reflected by an appropriate number of 
times from the basic segment, thus will coincide in ID and location with a unique instrumented 
point outside of the basic segment on the test article. A cloned point is a grid point in the basic 
segment with a unique grid point ID number, which has duplicate coordinates of a multiply- 
reflected direct or image point, but has no structural function. There are as many clones at a 
point as there are multiplicities of instrumented reflections to a point. Clones are fully constrain- 
ed points that have no elements connected to them. The point whose coordinates are being 
cloned is called a parent point. It is usual for the analysis to have more data points than the test 
data has, so it is necessary to partition the analytical mode shapes down to contain data corre- 
sponding to only those points that were instrumented on the test article. Partitioning vector 
TESET earmarks all direct, image, and cloned points so that when TESET acts on the matrix of 
analytical mode shapes, PHIG, the modal deflections for instrumented points only will be retained 
in the matrix PHITE. But PHITE, at this raw stage, is not ready for comparison with test data, be- 
cause modal data from the basic segment has to be reflected out to the instrumented points with 
correct signs. The signs of the deflections of instrumented points to be assigned to image and 
cloned points, depend on their location and the type of symmetry being simulated by an analyti- 
cal run. A way has been devised to communicate these particulars to the processor program which 
will enable it to convert PHITE to make it ready for STATCORR. The same parameters that de- 
scribed the type of symmetry for the case of both sources using symmetry will serve equally well 
for this case. Parameters SYM1PLAN, SYM2PLAN and SYM3PLAN must be supplied so that 
they will be output to the postprocessor. In order to so communicate a location, a convention 
must be adapted. The following is the convention used for identifying the octants. Octant # 1 
is the basic segment. Octant # 2 is the reflection of the basic segment about the first plane of sym- 
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metry. By-pass, momentarily, the definition of octant # 3 and define octant # 4 as the reflection 
of the basic segment about the second plane of symmetry. Octant #3 lies in between octants 2 and 
4. It is a double reflection of the basic segment about the first and second planes of symmetry. 
Octant # 5 is the reflection of the basic segment about the third plane of symmetry. Rather than 
define octants 6, 7, and 8 as a series of precise double and triple reflections, it is easier to picture 
them as lying under octants 2, 3, and 4, just as octant # 5 lies under the basic octant. It remains 
to set a procedure for reporting the locations of instrumented points. This is done by reorganiz- 
ing the data, already contained in the TESET vector, into an 8 columned matrix which will be 
named DOFLIST. Columns one through eight of DOFLIST refer to octant locations one through 
eight. Direct points are tabulated in ascending sequence of internal DOF number in column one. 
Image points and parents of cloned points reflected from octant # 2 are tabulated in ascending se- 
quence of internal DOF number in column two. Similarly, tabulations follow for octants 3 
through 8. Tabulations are omitted for any octant for which symmetry is absent (i.e., SYM#PLAN 
is zero). The essential difference between TE SET and DOFLIST will be stressed at this juncture. 
TESET is a vector tabulation of internal degrees of freedom vs. instrumented grid point numbers 
including clone IDs. DOFLIST is an 8 columned matrix tabulation of internal degrees of freedom 
vs. all instrumented grid point numbers except that parent IDs are used for clones. 

By definition, if SYM2PLAN is equal to —1, the reflection with respect to the second plane 
is anti-symmetric, and the vertical displacements in octant # 4, as shown in Figure 2, are of the op- 
posite sign from those in the basic segment. Similarly, if SYM3PLAN is equal to +1, the vertical 
displacements in octant 5 would also be of opposite sign from those of the basic segment. To de- 
termine the sign for octants # 3, 6, 7, and 8, the parameters must be used in particular ways. 
Thus, the sign of displacements in octant # 3 depend on the types of symmetry across both the 
first and second planes of symmetry. If one were symmetric and the other anti-symmetric, the 
vertical displacements for octant # 3 would be of opposite sign to the basic. This discussion deals 
with a given point in the basic segment which is paired with its image point in the reflected oc- 
tant; therefore all such pairs of displacements are equal in magnitude. Symmetry and anti-sym- 
metry for 3 translational and 3 rotational degrees of freedom will be defined with respect to a 
plane of symmetry, and then is translated to items of action pertaining to the correlation coeffici- 
ent under the option of “symmetry for analysis and whole structure for test.” Consider, for 
discussion purposes, that the plane of symmetry is a coordinate plane and, it is oriented vertically. 
Designate one of the other two coordinate planes to be the horizontal plane, and the remaining 
one to be transverse. 

For the moment, language will be adapted" for describing the modal displacements in terms 
of this convention for the 3 planes. Displacements perpendicular to the plane of symmetry can be 
described as being towards or away from the plane of symmetry. Displacements perpendicular 
to the horizontal plane can be described as being up or down. Displacements perpendicular to the 
transverse plane can be described as back and forth. The goal is to define the behavior of the re- 
flection with respect to the plane of symmetry for each of the 6 displacement components at a 
Grid Point, for either the symmetric or the anti-symmetric case. 

When data is collected on a whole structure with respect to a single fixed coordinate system 
with coordinate planes imbedded in the frame of the planes of symmetry, it is unaware of the de- 
scriptions: twd/away, up/down, back/forth. Therefore, when analytical symmetric data is convert- 
ed from the basic segment to modal behavior in the seven other octants, it must take into account 
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the way that displacements will appear in the fixed coordinate system of the whole. As an ex- 
ample, according to Table S, the symmetric behavior involving displacements perpendicular to the 
plane of symmetry will result in the basic and image points both moving synchronously toward or 
synchronously away from the plane of symmetry. In terms of the coordinate system of the 
whole, the displacements involve the same coordinate axis. Since they are both moving toward 
the plane of symmetry, they are moving towards each other; that is, they bear opposite signs. 


Table S 

IMAGE VS BASIC MOTIONS WITH RESPECT TO PLANES OF SYMMETRY 


DIRECTION 
OF MOTION 

TRANS i P/S 


ROT’Nl P/S 


RELATIONSHIPS 


Both twd v. „ 

„ > - Symmetric 

Both away 

One twd & one away - AntiSym 


Avert/ A forth: 


Both A up 
or Both Adown 
with Both A forth 


— Symmetric 


One Aup & one Adownv. . 
with Both A forth >-AntiSym 


TRANS 1 H/P 


ROT’N 1 H/P 


Both up v. c 

Both down> - Symn,etnc 

One up & one down — AntiSym 


Aacross/ A forth: 


Both A twd v c . 
Both 4away > ~ Symmelnc 
with Both A forth 
One A twd & one A away 
with Both A forth 


\ - AntiSym 


TRANS 1 T/P 


ROT’Nl T/P 


Both forth\ c . . 

Both back >- Symmetnc 

One forth & one back — AntiSym 

Aacross/Aup: Both Atwd 

or Both Aaway 
with Both Aup 
One Atwd & one Aaway 
with Both Aup 


— Symmetric 


^ AntiSym 


Consequently, adaptation of basic segment data to octant # 2 symmetrical modal behavior, re- 
quires that the sign be changed for those components which are perpendicular to the plane of sym- 
metry. Currently, it is recognized that only translational data is measured in test, but Table U will 
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be prepared for sign changes not only for the 3 translational components, but also for the three 
rotational components. In order to implement this conversion, certain restrictions must be placed 
on the analytical model. Only the basic coordinate system will be allowed to be used for grid 
points in the basic segment which correspond to instrumented test points. Its coordinate planes 
must be parallel to planes of reflective symmetry. To report how the basic coordinate system is 
oriented to the planes of symmetry, a vector is required to be constructed on DMI cards called 
IDCORD. The basic coordinate system ID (non-zero) is entered as an integer in the row position, 
and the translational component which is perpendicular to the lsJ plane of symmetry is entered 
as a real number, followed by the translational component, which is perpendicular to the 2nd 
plane of symmetry. The size of the matrix to be entered on the “0” card is: rows = ID # + 1, 
columns =1. If a small number is used for the coordinate system ID, it will eliminate the unneces- 
sary processing of zeroes. The card will have the appearance of that shown in Figure 3, page 18. 

Table U 


OCTANT 
LOCATION 
OF IMAGE 
POINT 


ACTION TABLE FOR CONVERTING BASIC SEGMENT DATA 
SIGN OF DISPLACEMENT 


TYPF. & SEQUENCE 

1 

2 3 4 

5 6 

OF REFLECTION 

2 

3 1 5 

6 4 

SYM 1 

SYM2 

SYM3 

3 

1 2 6 

4 5 

+ 1 



CHNG 

NOOP NOOP NOOP 

CHNG CHNG 

-1 



NOOP 

CHNG CHNG CHNG 

NOOP NOOP 


+ 1 
1 


NOOP 

CHNG NOOP CHNG 

NOOP CHNG 



CHNC, NOOP CHNC. NOOP 

CHNG NOOP 



-Hi 

NOOP 

NOOP CHNG CHNG CHNG NOOP 



-1 

CHNG 

CHNC. NOOP NOOP 

NOOP CHNG 

+ 1 
1 

1 


CHNG 

NOOP NOOP 

CHNG 

+ 1 


NOOP 

CHNG CHNG NOOP 

+ 1 
I 

+ 1 



NOOP 

CHNG 

- 1 



CHNG 

NOOP 

+ 1 
1 


- 1 

CHNG 

NOOP NOOP 

CHNG 


+ 1 

NOOP 

CHNG CHNG 

NOOP 

+ 1 
1 


+ 1 


NOOP 

CHNG 


-1 


CHNG 

NOOP 


+ 1 
I 

- 1 


CHNG NOOP 

NOOP CHNG 


+ 1 


NOOP CHNG 

CHNG NOOP 


+ 1 

+ 1 

NOOP 

CHNG 



- 1 

1 

CHNC, 

NOOP 


+ 1 
1 

I 

1 

CHNG 

NOOP 


+ 1 
1 

+ 1 

NOOP 

CHNG 


+ 1 
1 

+ 1 


NOOP 

CHNG 

+ 1 

1 


CHNG 

NOOP 

+ 1 

+ 1 
1 

1 


NOOP 

CHNG 

1 

+ 1 


CHNG 

NOOP 


COMP PERP TO TP P/S 
1ST 2ND 3RD 


3 
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Notice that the octants are not entered in ascending sequence. The scheme used here was to work 
towards increasing complexity starting with octants involving only a single reflection, then double, 
and Finally triple reflection. 


DMI IDCORD 0 

DMI IDCORD 1 


2 

ID 

(INTEGER) 


1 

PERP 1 ST 
(REAL) 


1 

PERP 2ND 
(REAL) 


ID + 1 1 


Figure 3 


What action needs to be taken to convert a component in the basic segment to a modal displace- 
ment for an image point in an external octant has been worked out for either type of symmetry, 
for all six component displacements, and for the three distinct orientations of the basic coordinate 
system. The action to be taken is, either to change the sign of that component or not. If the sign 
is to be changed, the notation in the table is CHNG. If the sign is to be left unchanged, the entry 
in the table is expressed in computer jargon as NOOP, meaning no operation. 

Certain components can participate in only certain types of multiple reflections. Null entries 
are mutually exclusive. The processor program will accomplish the actions tabulated here, so that 
the analytical mode shapes delivered to the STATCORR program will represent whole structure 
behavior at those image points corresponding to instrumented points outside the basic segment. 


CODES 

DMAP ALTER’s; (STATDMAP) 

In reviewing the design of the DMAP ALTER packet, one sees that its principal function is 
to transfer important data blocks from the NASTRAN runs to the STATCORR program. It uses 
the punch File for all output. It manages the partitioning of the eigenvectors to the order of test 
modes. 

A DMAP ALTER packet must be included with each analytical run in order to provide the 
data sets and parameters necessary to prepare the output properly for operating the STATCORR 
program. A listing of the ALTER packet is contained in this report on pages 22 through 26. 
As of this writing, the packet is compatible with NASTRAN Rigid Format Series P. It was intended 
to be self explanatory by supplementing it with comment statements. It is stored in the user library 
account ULIB1 in the Design Engineering Branch under the filename DRAO: [NASCAT. UTILITY] 
STATDMAP.LIS;! . Analysts who wish to implement this ALTER can access this File, then copy 
and merge it into the executive control of their own jobs. Scanning the left hand margins of 
STATDMAP, the reader will notice that there are several subdivisions. Each subdivision will be 
called a packet according to an identifier. Each ALTER packet subdividion is readily distinguished 
by a title. The first packet is entitled FREQUENCY DATA and uses the margin identifier $F. 
PARTITIONED MODE SHAPES is the title of the second packet, identified by $P. The third 
packet is entitled SYMMETRY FOR BOTH, using $S at the margin. Finally, the ANALYTICAL 
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SYMMETRY ONLY packet uses identifier $W. Packets F and P are always used. When data 
comes from a basic segment of symmetry for both analysis and test, packets F, P, and S are used. 
The entire set of ALTER packets is needed when the test data is based on the whole structure, but 
analysis uses a basic segment only. Packet $F delivers the LAMA table from a NASTRAN eigen- 
value execution to the punched output file using the DMAP utility module, TABPCH. Although 
the entire LAMA table is exported, the preprocessor program selects only three entries for each 
mode. These are words 5, 6, and 7 in Record 2 of LAMA giving the cyclic frequency, generalized 
mass, and generalized stiffness, respectively. Packet $P instructs the user how to set up his DMI 
cards for defining the partitioning vector TESET. Module PARTN employs TESET to partition 
the matrix of eigenvalues PHIG from the analytical size down to experimental size. Module OUT- 
PUT 3 delivers test-sized eigenvectors PHITE and vector TESET to the punch output file in DMI 
format. 

Packet $S provides a place for the analyst to catalog which of the planes of symmetry he 
named first, second, and third. There are options for rectangular, cylindrical or spherical coordi- 
nate systems, so the user is advised to erase the unused options to avoid the clutter. The primary 
purpose of packet $S is to define the type of reflective symmetry which was simulated for a given 
computer run. Each time a run is made with a different set of boundary conditions at the planes 
of symmetry, the appropriate entries should be made for the SYM # PLAN parameters. These 
parameters are not used at all in the NASTRAN execution. DMAP is being used as a transfer agent 
to pass the values of these parameters into the post-processor program in a self-contained package 
so that none of the pertinent data goes adrift or gets mixed with the wrong set. Module PARAM is 
set up for the user to input their values, and module PRTPARM transfers them to the output file. 
All other previous output has been sent to the punch output file, and these parameters should be 
no exception. This is done by intercepting the print output from PRTPARM just long enough to 
route it to the FORTRAN logical unit for punch output, and then restoring it. PARAM statement 
S4 directs the parameters to the punch output, and S7 restores it to print output. 

Packet $W instructs the user how to set up his DMI cards for tabulating the image points in 
the octants for matrix DOFLIST and for the IDCORD card. Module OUTPUT3 delivers the 
DOFLIST matrix and the IDCORD vector to the punch output file in DMI format. Once again, 
DMAP is being used strictly as a transfer agent to pass the data from the user to the processor pro- 
gram. 

Statement FI is inserted at position 103 of RF3 after the READ module has output LAMA. 
It picks up data block LAMA, which is a table, and sends it to the punch output file in DTI for- 
mat. Parameter 1 is given the BCD value “FN” to use in setting up the continuation cards. 

Statement PI does a row partition on the matrix of eigenvectors, PHIG. It is inserted at 
R.F. 3 statement #111 after the matrix PHIG has been checkpointed. The position of the col- 
umn partitioning vector is left deliberately unspecified, while TESET is entered at the position of 
the row partitioning vector (RP). TESET has values at internal sequence positions of the G set 
which correspond to the degrees of freedom of points that were instrumented on the test article. 
PARTN separates the rows into two sub matrices according to whether there is a value in the par- 
titioning vector for that row or not. The symbolic name of the partition matching the row of RP 
without values is Aj j , and that, matching the rows of RP with values is A 21 . Since our interest is 
in the A 21 partition only, all positions of the output data blocks from PARTN are left purged, ex- 
cept the second. The name given to the A 21 partition is PHITE. Giving a positive value of the 
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parameter SYM in the first position indicates a non-symmetric partition. The positive value given 
to parameter, TYPE, in the second position indicates that the output partitions are real single pre- 
cision matrices. PHIG is a real single precision matrix so any of its partitions will be the same. The 
value 2 for the 3rd and 4th positions indicates that their forms will be general rectangular. 

Statement P2 picks up data blocks PHITE and TESET and writes them on the punch output 
file in DMI format. Parameter PI in the first position is negative to indicate that data being output 
should be sent to the FORTRAN logical unit. BCD value, MOD, is assigned to the second parame- 
ter for use in setting up continuation cards for the PHITE matrix, and BCD value, DOF, is assigned 
to the third parameter for use in setting up continuation cards for the TESET vector. 

The explanation for why ALTER 125 was selected for inserting packet P will be delayed un- 
til after the explanation of statement S8. Statements SI, S2, and S3 define parameters SYM#PLAN 
through the module PARAM using logical operator MPY (multiply). In each instance, variable 
parameter INI is multiplied by the constant 1 to produce the output value. In effect, the value 
given to the variable parameter INI is the value of the output. 

Statement S4 uses the module PARAM to change the value of the parameter which controls 
the printing of NASTRAN output (OUTTAP) that is located in the system common block (SYS- 
TEM) of NASTRAN. The purpose of S4 is to change the OUTTAP assignment from the value of 
the FORTRAN logical unit for printing to the value of the FORTRAN logical unit for punching. 
The logical operator SYST uses the value for the 3rd parameter to catalog the sequence number of 
of the parameter in the system common block. Referring to Section 2. 4. 1.8 of the Programmer’s 
Manual, one Finds that the sequence number for parameter OUTTAP is 2, so 2 is entered as the 
value for the third parameter. SYST uses the value of the 4th parameter to catalog the assignment 
number of the FORTRAN logical unit. These FORTRAN assignments are arbitrary and depend 
upon the decision of the systems programmer at the time that NASTRAN is installed on a com- 
puter. The punch output was assigned to FORTRAN logical unit 77 at the time that NASTRAN 
was installed on the computer used for developing this method, so in order to direct the output 
data blocks of a module under the management of system common block, named OUTTAP, to go 
to the punch buffer, 77 is entered as the value for the 4th parameter. The BCD value DUM given 
to the 2nd parameter is meaningless, because this parameter is not used in the operation of the 
logical operator SYST. 

Statements S5, S6, and S7 involve module PRTPARM. The value 0 assigned to the 1st 
parameter enables the module to operate as a parameter printer. The particular BCD values as- 
signed to the 2nd parameter discriminate amongst all possible parameters and will select only 
those having the indicated BCD names. 

Statement S8 uses the module PARAM in the same sense as in statement S4 to assign the 
FORTRAN logical unit 6 to the system common block parameter OUTTAP. Making this tem- 
porary reassignment of a major regulation item should be treated gingerly so as not to interfere 
with normal operations. Rapid restorationn is imperative. Care must be exercised when PARAM 
logical operator SYST is put into operation. Since the printer output block is involved, care must 
be taken to ensure that anything destined for the print buffer has been already sent there, so that; 
only those things that are governed by parameter OUTTAP during this temporary assignment are 
routed to the punch buffer. Statement 125 occurs after module OFP has directed OEIGS, LAMA, 
and DQM1 to the print buffer and before new quantities from SDR2 are being readied for print- 
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ing. Thus diverting OUTTAP temporarily to punch at statement 125 does not endanger the print- 
ing of any important quantities. 

Statement W1 is inserted at position 128 to isolate it from those in packet S, as well as, to 
make sure that it would not be caught in a bypass caused by a conditional jump. It picks up the 
data blocks delivered by the analyst in the bulk data in DMI format, and writes them on the 
punch output file in DMI format once again. Parameter PI in the 1st position is negative to in- 
dicate that the data output should be sent to the FORTRAN logical unit. The BCD value, LIS, 
is assigned to the 2nd parameter, and FRM to the 3rd parameter for use in setting up the continu- 
ation cards. 
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• STATCORR PREP 

t DMAP ALTER PACKET 

S 

• 

9 FOR USE WITH RIGID FORMAT SERIES P 

9 

9 

9 

9 

ALTER 103 

• THIS DMAP ALTER PACKET IS DESICNED TO PROVIDE OUTPUT TO PROCESSORS 

• WHICH IN TURK PREPARE THE ANALYTICAL INPUT TO THE STATISTICAL COR- 

• RELATION PROCRAM ***STATCOIUl*** FOR COMPARING ANALYTICALLY DERIVED 

• MODAL DATA WITH EXPERIMENTALLY DERIVED MODAL DATA. 

• SUBPACKETS ***F*** AND ***P*** WILL ALWAYS BE USED. SUBPACKET ***S*** 
f IS DESIGNED FOR USE WHEN BOTH ANALYTICAL AND EXPERIMENTAL DATA ARE 

t REPORTED FOR ONLY A SINCLE SECMENT OF REFLECTIVE SYMMETRY IN ONE, 
t TWO, OR THREE COORDINATE DIRECTIONS. SUBPACKET ***W***IS DESIGNED 

• FOR THE CASE WHEN REFLECTIVE SYMMETRY IS USED WITH .ANALYSIS BUT NOT 

• WITH TEST. 

9 

9 

9 

f*******************FREQUENCY DATA*********************************** 

9 

•FI 

TABPCH LAMA, , , ,/VFN • 

• 

•F WORD 5 IN EACH ENTRY OF RECORD 2 OF THE LAMA TABLE GIVES THE CYCLIC 
•F FREQUENCY FOR A MODE. 

•F WORD 6 IN EACH ENTRY OF RECORD 2 OF THE LAMA TABLE GIVES THE GENER- 
•F ALIZED MASS FOR A MODE. 

•F WORD 7 IN EACH ENTRY OF RECORD 2 OF THE LAMA TABLE CIVES THE GENER- 
•F ALIZED STIFFNESS FOR A MODE. 

• 

•*****************PARTITIONED MODE SHAPES*************************** 

•P 

ALTER 1 1 1 

••P 

•P PHIC IS PARTITIONED TO THAT SET OF FREEDOMS WHICH ARE TO BE 
•P INSTRUMENTED ON THE TEST ARTICLE. 

•P 

•P THE PARTTIONINC VECTOR TESET IS INPUT USING DMI . THE PRECEDING 
•P EIGENVALUE RUN SHOULD INCLUDE DIAG 21,22 TO GIVE THE CORRELATION 
•P BETWEEN THE INTERNAL D.O.F.'S AND THE EXTERNAL NUMBERING. THE 
•P PARTITIONING VECTOR IS ORGANIZED ACCORDING TO INTERNAL SEQUENCING. 

•P MODULE PARTN IS ORGANIZED TO DO A ROW PARTITION AMONGST THE DEGREES 
•P OF FREEDOM IN ORDER TO PRESENT A DESCRIPTION OF THE MODAL BEHAVIOR 
•P ONLY IN TERMS OF THE INSTRUMENTED D.O.F.’S. 

•P 

•P THE ZERO CARD OF DMI HAS THESE ENTRIES: FIELD 1, DMI » FIELD 2, TESET » 
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•P FIELD 3 , ZERO 

•P"V I ELD 4, 2 — TO INDICATE THAT THE MATRIX BEING BUILT IS RECTANGULAR; 
•P^riELDS 5 8 6, 1— TO INDICATE SINGLE PRECISION IN AND SINGLE PRECISION 
•P OUT; FIELD 7, BLANK; FIELD 8, G-SIZE; FIELD 9, 1— TO INDICATE A 
S> SINGLE COLUMN OUTPUT. 

•P $P P 

•P THE NEXT DATA CARD OF DMI HAS THESE ENTRIES: FIELD 1, DM I ; FIELD 2, 

•P TESET; FIELD 3, 1 — TO INDICATE DATA FOLLOWING PERTAINS TO COLUMN 1 
•P OF THE OUTPUT MATRIX; 

•P FIELD 4, THE INTERNAL SEQUENCE NUMBER, IN INTEGER FORMAT, OF 
•P THE FIRST INSTRUMENTED D.O.F. ACCORDING TO THE DIAC 21 TABULATION; 

•P FIELD 5, THE EXTERNAL GRID POINT ID NUMBER CORRESPONDING TO THE 
•P FIRST DOF. IT IS TO BE WRITTEN AS A REAL NUMBER, I.E. GP ID 
fP FOLLOWED BY A PERIOD. 

•P FIELD 6, THE SECOND INSTRUMENTED D.O.F. PER DIAG 21; INTEGER. 

•P FIELD 7, THE EXTERNAL GRID POINT ID NUMBER CORRESPONDING TO THE 
•P SECOND DOF. IT IS TO BE WRITTEN AS A REAL NUMBER. 

•P CONTINUE IN A SIMILAR MANNER UNTIL ALL INSTRUMENTED POINTS HAVE 
•P BEEN ENTERED. IN THE CASE OF SUBPACKET ***W*** THE GRID POINT 
•P NUMBERS INCLUDE CLONES AND IMAGES. 

•PI 

PARTN PHIG, , TESET/, PHITE, ,/C,N,+l/C,N, 1/C, N, 2/C, N, 2 S 

•P 

•P MATRIX PIIITE IIAS BEEN PARTITIONED FROM THE PIIIG DATA BLOCK. PIIITE 
•P CONTAINS MODAL DATA FOR ONLY THOSE DOF THAT WERE INSTRUMENTED. 8 

•P 

•P2 

0UTPUT3 PHITE, TESET //C, N, - 1/C, Y, Nl= VCT /C,Y,N2=D0F • 

•P 

•P PHITE AND TESET ARE OUTPUT IN DMI FORMAT AS PUNCH FILES. 

•P IF THE USER WANTS THE 

•P DMI CARD IMAGES TO BE WRITTEN TO A DISK FILE INSTEAD OF BEING PUNCHED, 
•P THE USER MUST SUPPLY THE APPROPRIATE JOB CONTROL LANGUAGE STATEMENTS. 

• 

#**************************symmetry for both*************************** 

• 

9 

9S THE ANALYST SHOULD SUBSTITUTE THE PROPER ENTRIES FOR EACH OF THE 
•S VARIABLES REPRESENTED IN *X*X*X*X* FORMAT JN THE NEXT PARACllAPU. 

• 

•S 

•S IF THE ANALYSIS IS TO BE CONFINED TO ORLY ONE SEGMENT OF REFLECTIVE 
•S SYMMETRY AND SIMULTANEOUSLY IF TEST RESULTS WILL ALSO BE REPORTED 
•S FOR ONLY ONE SEGMENT, THEN THE TYPE OF MODES BEING PRODUCED IN EACH 
fS RUN MUST BE INDICATED IN THE FOLLOWING THREE PARAMETER STATEMENTS. 

•S FOR RECORD PURPOSES, THE FIRST PLANE OF SYMMETRY IS PERPENDICULAR 
•S TO THE *A*A*A*A* AXIS AND CUTS THAT AXIS AT *A*A*A*A* (RECTANGULAR), 
tS OR IS AN R: Z PLANE THRU THETA = *T*T*T*T* , PHI = 0 (CYLINDRICAL), 

•S OR IS A PLANE THRU THETA = 6 , PHI = 0 (SPHERICAL). 
tS THE SECOND PLANE OF SYMMETRY IS PERPENDICULAR TO TIIE *B*B*B*B* AXIS 
•S AND CUTS THAT AXIS AT *B*B*B*B* (RECTANGULAR), OR IS AN R:Z PLANE 
tS THRU THETA = [*T*T*T*T* + PI/2] , PHI = 0( CYLINDRICAL) , OR IS A PLANE 
•S THRU THETA = PI/2, PHI = 0 (SPHERICAL). 

•S THE THIRD PLANE OF SYMMETRY IS PERPENDICULAR TO THE *C*C*C*C* AXIS 
•S AND CUTS THAT AXIS AT *C*C*C*C* (RECTANGULAR), OR IS PERPENDICULAR 
•S TO THE Z AXIS AND CUTS THAT AXIS AT Z = *V*V*V*V* (CYLINDRICAL), 

•S OR IS A PLANE AT PHI = PI/2 FOR ALL THETA (SPHERICAL). 
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ss 

9 

9S 

SS THE DESIGNATION FOR A SYMMETRICAL REFLECTION W.R.T. AN AXIS IS +1; 

SS ENTER INI = +1 ON THE PARAII CARD IN PLACE OF XXX. 

SS THE DESIGNATION FOR ANT I -SYMMETRIC REFLECTION W.R.T. AN AXIS IS -1; 

SS ENTER INI = -1 ON THE PARAM CARD IN PLACE OF XXX. 

SS THE DESIGNATION FOR ABSENCE OF SYMMETRY WITH RESPECT TO AN AXIS IS 0; 
SS ENTER INI = 0 ON THE PARAM CARD IN PLACE OF XXX. 

$ ~ 

SS 

SSSS ALTER !25 
SSI 
S 
9 


SS VALUES CV SYM'PLAlf WILL CHANCE EVERY TIME THE CONSTRAINTS 
SS OF SYMMETRY CHANGE. 

S 

S 

S 

SS1SPARAM //C, N, MPY/V, Y, SYM1PLAN/V, N, IN1=XXX/V,N, IN2= 1 3 

SS 

SS THIS STATEMENT SAYS WHETHER THE MODE SHAPE IS SYMMETRIC OR ANTI- 
SS SYMMETRIC OR HAS NO SYMMETRY WITH RESPECT TO PLANE I. 

C3 

552 

SS2GPARAM //C, N, MPY/V, Y, SYM2PLAN/V, N, IN1 = XXX/V,N, IN2= 1 9 
SS 

SS THIS STATEMENT TELLS WHAT KIND OF SYMMETRY PERTAINS TO PLANE 2. 

SS 

553 

SS33PARAH //C, N , MPY/V, Y, SYM3PLAN/V, H, IN1 S XXX/V, N, IN2= 1 3 
SS 

SS THIS STATEMENT TELLS WHAT KIND OF SYMMETRY PERTAINS TO PLANE 3. 

S 

C PARAII //C,N,SYST/V,N,DUM/C,N,2/C,W,77 9 
SS 

SS CnANCES THE FORTRAN LOCI CAL UNIT FOR TIIE SYSTEM COMMON DATA BLOCK, 
SS CALLED OUTPUT, FROM ITS USUAL VAX VALUE 

SS OF 6 TO '• i.E VAX VALUE RESERVED FOR PUNCH WHICH IS 77, SO THAT IT CAN 
SS BE WRITTEN TO THE SAME DISK FILE AS THE CHKPNT DICTIONARY AND THE 
SS 0UTPUT3 FILES. 

SS 


555 

SS5GPRTPARM //C f H f 0/C, N. SYM1 PLAN 9 
SS 

556 

SS6SPRTPAIU I //C, N , 0/C , N , S YM2PLAN 9 

SS 

557 

SS7SPRTPARM //C , N , 0/C , N t S YM3PLAN 9 

SS 

ssa 

SSHSPARAM //C, N, SYST/V, N, RUM/C, N, 2/C, N, 6 9 
SS 

$S RESTORED HIE FORTRAN LOGICAL UNIT FOR OUTPUT BACK TO THE VAX VALUE 
SS OF 6 NOW THAT THE NEED FOR ASSIGNING IT TO 77 IS OVER. 

8 

9 
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• 

S** fc**************ANALYTICAL SYMMETRY ONLYfcsMe*********************** 

SW 

SW JSE SUBPACKET **W** IF EXPERIMENTAL DATA COMES FROM TIIE WHOLE STRUCTURE 

SW BUT ANALYTICAL DATA COMES FROM A BASIC SEGMENT OF SYMMETRY . TIIE 

SW ANALYTICAL BASIC SEGMENT WILL CONTAIN POINTS CORRESPONDING TO TUOSE 

SW TEST POINTS IN THE BASIC SEGMENT. IN ADDITION THE ANALYTICAL nASIC 

SW SEGMENT WILL CONTAIN POINTS WHICH ARE REFLECTED IMAGES OF ALL TEST POINTS 

SW WHICH ARE LOCATED OUTS IDE OF THE BASIC SEGMENT. THE PARTITIONING 

•W VECTOR ***TESET*** WILL INCLUDE THE DIRECT POINTS PLUS ALL CLONED AND 

SW IMAGE POINTS. ALL REFLECTED PO INI'S 

•W MUST BE TALLIED FOR ***STATCORR*** TO PROCESS THE DATA 

•W CORRECTLY. THIS TALLY IS COMMUNICATED IN THE FORM OF AN EICIIT 

•W COLUMNED MATRIX CALLED ***DOFLIST*** ■*- 

SW 

SW 

SW 

SW 

SW WARNING. WARNING. BE SURE THAT CLONED CRID POINTS 

SW ARE INCLUDED IN THE BULK DATA AND ARE FULLY CONSTRAINED. 

SW 

SW 

SW 

SW 

SW EACH OF THE 8 COLUMNS IS RESERVED FOR A PARTICULAR OCTANT. THESE 
SW ARE THE DEFINITIONS OF THE OCTANTS. 

SW OCTANT 1 IS THE BASIC SECMENT. OCTANT TWO IS THE REFLECTION OF THE 
SW BASIC SEGMENT ABOUT SYM1PLAN. OCTANT FOUR IS THE REFLECTION OF THE 
SW BASIC SEGMENT ABOUT SYM2PLAN. OCTANT THREE IS IN BETWEEN OCTANTS 
SW TWO AND FOUR. OCTANT FIVE IS THE REFLECTION OF TIIE BASIC SEGMENT 
SW ABOUT SYM3PLAN. OCTANTS SIX, SEVEN, AND EIGHT ARE ROTATED AWAY 
SW FROM OCTANT FIVE IN THE SAME SENSE THAT OCTANTS TWO, THREE, AND 
SW FOUR ARE ROTATED AWAY FROM OCTANT ONE. 

SW DATA FOR TALLYING THE IMAGE POINTS REFLECTED FROM EACH OCTANT IS 
SW ORGANIZED IN DM! INPUT. COLUMN ONE IS RESERVED FOR LISTING ALL 
SW POINTS PLUS COMPONENT DIRECTIONS WHICH ARE COMPANIONS 
SW OF THE INSTRUMENTED POINTS IN THE 

SW BASIC SEGMENT. COLUMN TWO IS RESERVED FOR LISTING ALL UNIQUE 
SW POINTS OR PARENTS OF CLONES 

SW IN THE BASIC SECMENT WHICH ARE REFLECTIONS OF INSTRUMENTED 
SW POINTS FROM OCTANT * 2 . SIMILARLY, COLUMN THREE IS RESERVED FOR 
SW LISTINC ALL UNIQUE POINTS OR PARENTS OF CLONES IN THE BASIC SEGMENT 
SW WHICH ARE REFLECTIONS OF 

SW INSTRUMENTED POINTS FROM OCTANT * 3 . COL 4 LIST FOR OCTANT 4 AND 
SW SO FORTH THRU COL 8 FOR OCTANT 8. 

SW EXAMPLES OF EACH OF THE NINE DMI LOGICAL CARDS FOR THIS MATRIX 
SW EXPRESSED IN FREE FIELD FORMAT ARE: 

SW DMI ,DOFLIST,0,2, 1, 1 , BLANK, G-S IZE, 8, BLANK 

SW DMI , DOFLIST, 1 , INT DOF *,GP *, INT DOF *,GP *, INT DOF *,CP *,CONT’N 

SW DMI .DOFLIST, 2, INT DOF #,GP *, INT DOF #,GP *, INT DOF *,GP *,CONT’N 

SW DMI DOFLIST, 3, INT DOF *,GP *, INT DOF *,GP * , INT DOF *,GP *,CONT’N 

SW DMI ETC FOR COLUMNS 4 THRU 8. 

SW NOTICE THAT THE CONTENT OF CARDS FROM COLUMN ONE ET SEQ. AnE 

SW PAIRED QUANTITIES. THE FIRST OF THE PAIR IS THE ROW NUMBER 

SW WHICH IS THE INTERNAL DECREE OF FREEDOM NUMBER WRITTEN AS AN 
SW INTEGER. THE SECOND OF THE PAIR IS A REAL NUMBER FORM OF 
SW THE GRID POINT ID CORRESPONDING TO THE INTERNAL DOF . 

SW PAIRS FOR A GIVEN COLUMN 

SW OVERFLOW ONTO CONTINUATION CARDS UNTIL ALL POINTS FOR AN OCTANT 
SW ARE TALLIED. DATA FOR THE SUCCEEDING COLUMN CONSTITUTE A NEW 
SW LOGICAL CARD STARTING WITH THE COLUMN NUMBER THEN FOLLOWED BY ITS 
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*V SUBSEQUENT PAIRS. 

tV .... * • , 

tv 

tv 

tv 

tV A SECOND Dm VEC TOR VA RIED *:mDCORD*** IS PUT INTO THE BEK DATA 
tV TO SPECIFY THE ORIENTATION OF BASIC COORDINATE P LAN ES VI TH TOE 

tv planes of syHMTKy. " "ValM #T t& far bas ic c ompon ent that rs 

tv PERPENDICULAR TO TOE 1ST PLANE OF SYMMETRY. VALUE *2 IS TOE 
tV BASIC COMPONENT THAT IS PERPENDICULAR TO THE 2ND PLANE OF 
tV SYMMETRY. SET VALE *2 = 0.0 IF THERE IS ONLY ONE PLANE OF 
tV SYMMETRY. NO ENTRY IS MADE FOR THE THIRD COMPONENT REGARDLESS 
tV OF TOE EXISTENCE OF A THIRD PLANE OF SYMMETRY. 
tV 
tV 

tv 

tv VS ALTER 126 
tVl 

tVlSOUTPUTO DOFLIST, IDCORD//C, N,- 1/C, Y, LIS/C, Y,FRM S 

tV POST PROCESSOR VILL APPLY TOE APPROPRIATE SIGN TO THE 
tV REFLECTED DISPLACEMENTS 

tV ACCORDING TO THE COMBINATION OF INI PARAMETERS ENTERED IN STATEMENTS 
$W SI, S2 , AND S3 . 

t 

END ALTER 


MERGE 


In case eigenvectors have been obtained from several disjoint NASTRAN runs, it is possible 
to combine these vectors into a single set for use with STATCORR. Probably some overlap 
of frequency ranges would occur with consequent duplication of vectors. The pre-processor 
program MERGE accepts Files of eigenvectors from several runs and discards all duplicates. It 
reassembles the resulting into a sequence according to ascending frequencies and delivers them 
as a single file. MERGE operates interactively and the user responds to the on-line prompts. 
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100 c 

200 G 
300 0 

400 C 
500 C 
600 C 
700 C 
800 G 
900 G 
1000 G 
1100 G 
1200 G 
1300 C 
1400 G 
1500 G 
1600 C 
1700 C 
1800 G 
1900 C 
2000 G 
2100 C 
2200 G 
2300 C 
2400 G 
2500 G 
2600 G 


PROGRAM MERGE 


THIS PROGRAM WILL MERGE TWO MATRIX FILES WITH THE SAME NUMBER 
OF COLUMNS INTO A SINGLE MATRIX WITH THE SAME NUMBER OF COL- 
UMNS AND A TOT AL NUMBER OF ROWS EQUAL TO THE SUM OF THE NUMBER 
OF ROWS IN THE ORIGINAL TWO MATRICES. 

A CORRESPONDING COLUMN IN EACH MATRIX (NOT NECESSARILY THE SAME 
COLUMN NUMBER) MU ST BE IN ASCENDING SORT, AND A SORTED MERGE WILL 
BE DONE BASED ON THE RESPECTIVE CORRESPONDING COLUMNS. 

IF DESIRED, DUPLICATED ROWS WILL BE DELETED; A MAXIMUM RELATIVE 
DIFFERENCE (EPSILON) MAY BE SPECIFIED FOR DEFINING "EQUALITY" ' 

OR A DEFAULT VALUE MAY BE USED. 

IF DESIRED, ROWS CONTAINING TERMS BELOW A MAXIMUM THRESHOLD VALUE 
WILL BE DELETED; THE THRESHOLD VALUE MAY BE SPECIFIED, OR A 
DEFAULT VALUE MAY BE USED. 

IF DESIRED, ROWS CONTAINING TERMS ABOVE A MINIMUM CEILING VALUE 
WILL BE DELETED; THE CEILING VALUE MAY BE SPECIFIED, OR A 
DEFAULT VALUE MAY BE USED. 

ADDITIONAL MATRIX PAIRS THAT ARE ASSOCIATED RESPECTIVELY WITH I THE 
ORIGINAL MATRIX PAI R MAY ALSO BE MERGED TOGETHER, BASED ON THE 
SORTED MERGE OF THE ORIGINAL PAIR. IN THIS CASE, COLUMNS OF AN 
ASSOCIATED MATRIX MUST CORRESPOND TO ROWS OF THE ORIGINAL MATRIX 
WITH WHICH IT IS ASSOCIATED. ANY DELETED DUPLICATE, BELOW- 
THRESHOLD, OR ABOVE-CEILING ROWS OF THE ORIGINAL MATRIX WILL 
RESULT IN CORRESPONDINGLY DELETED COLUMNS OF THE ASSOCIATED 
MERGED MATRIX. 
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PROGRAM MERGE 
C 

CHARACTER* 8 MATRIX1 , MATRIX2 
CHARACTER* 1 YESNO 
REAL SORTL I ST( 2 , 100) 

INTEGER MTXKEYC 200) , ROWKEY( 200) 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

100 

150 

C 

1000 

2000 


DO 100 ISET= 1 , 1000000 

ASK IF ANOTHER SET IS DESIRED 
TYPE 1000 

ACCEPT 2000, YESNO 
IF (YESNO .NE. 9 Y* ) GOTO 150 
REQUEST NAMES 0 OPEN MATRIX FILES (LOG, UNIT NOS. 102) 

CALL GETKEYMTX( MATRIX1 , MATRIX2, KROWl , KR0W2 , 

: KCOL 1 , KC0L2 , * 100) 

REQUEST KEY COLUMNS FOR SORTED MERGE, 0 READ THEM INTO ARRAY 
CALL GETSORTCOL( MATRIX1 , MATRIX2, KROWl , KROW2 , 

: KCOL 1 , KC0L2 1 SORTL I ST) 

GENERATE KEY LISTS FOR MERGING 

CALL KEYLIST( SORTL 1ST, KROWl , KR0W2 , MTXKEY, ROWKEY) 

SET TOTAL NUMBER OF ROWS (PRIOR TO DELETIONS) 

KOUT = KROWl + KR0W2 
CULL DUPLICATE VALUES 

CALL DELDUPES( SORTLIST, MTXKEY, ROWKEY, KROWl , KR0W2 , KOUTX 
CULL BELOW- THRESHOLD VALUES 

CALL DELTHRESH( SORTLIST, MTXKEY, ROWKEY, 
t KROWl , KR0W2 , KOUT, * 1 00) 

CULL ABOVE-CEILING VALUES 

CALL DELCEILING( SORTL 1ST, MTXKEY, ROWKEY, 

K KROWl , KR0W2 , KOUT ,*100) 

MERGE KEY MATRICES INTO OUTPUT MATRIX. 0 CLOSE MATRIX FILES 
CALL MERGEKEY( MTXKEY, KOUT) 

MERGE ASSOCIATED MATRIX PAIRS 

CALL ASSOC I ATE( MTXKEY, KROWl, KR0W2, KOUT) 

CONTINUE 

STOP ’END OF MATRIX SET REQUESTS’ 

FORMAT( / , * MERGE ANOTHER KEY MATRIX SET? ( Y OR N) : ’ ,•) 

FORMAT! Al) 

END 
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SUBROUTINE ASSOC I ATE ( MTXKEY, KROW1 , KR0W2, KDUT) 

C 

C MERGE ASSOCIATED MATRIX PAIRS 
C 

INTEGER MTXKEYC 200) 

CHARACTER* 1 YESNO 
C 

DO 100 IASSOC= 1 , 1000000 

C ASK IF ANOTHER ASSOCIATED MATRIX PAIR IS DESIRED 

TYPE 1000 

ACCEPT 2000, YESNO 

IF (YESNO .NE. ’ Y’ ) GOTO 150 

C OPEN ASSOCIATED MATRIX PAIR (LOGICAL UNIT NOS. 182) AND 

C CHECK FOR VALID DIMENSIONS. 

CALL GETASCMTX( KR0W1 , KR0N2, MROW, NCOL, * 100) 

C MERGE ASOCIATED PAIR INTO OUTPUT MATRIX FILE 8 CLOSE FILES 

CALL MERGEASC( MTXKEY, KOUT, MROW, NCOL) 

100 CONTINUE 

C NO MORE ASSOCIATED MATRIX REQUESTS 

150 TYPE *, * END OF ASSOCIATED MATRIX REQUESTS’ 

RETURN 

C 

1000 FORMAT!/, f MERGE ANOTHER ASSOCIATED PAIR OF MATRICES? (Y OR N) 

* , 0 ) 

2000 FORMAT! Al) 

END 
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SUBROUTINE DELDUPES( SORTLIST, MTXKEY, ROWKEY, KROW1 , KR0W2 1 KOUT) 
C 

C DELETE ROWS WHERE THE MERGING TERM IS A DUPLICATE (WITHIN 
C RELATIVE DIFFERENCE EPSILON) OF THE NEXT TERM. 

C 

INTEGER MTXKEYC 200) , ROWKEYC 200) 

REAL S0RTLISTC2, 100) 

CHARACTER* 1 YESNO 
C 

C ASK IF DUPLICATE DELETION IS DESIRED 
TYPE 1000 

ACCEPT 2000, YESNO 
IF (YESNO .NE. ’Y’) RETURN 

C SET DEFAULT EPSILON 8 REQUEST INPUT VALUE TO REPLACE IT 
EPS = 1.0E-5 
TYPE 3000, EPS 
ACCEPT *, EPS 

C CULL OUT DUPLICATE VALUES 
KOUT = KROW1 + KR0W2 
DO 1O0 1= 1,KR0W1+KR0W2-1 

VTHIS = SORTL IST< MTXKEYC I ) , ROWKEYC I ) ) 

VNEXT = SORTL I ST< MTXKEYC 1+1) , ROWKEYC I+D) 

DIFF = ABSC VNEXT- VTHIS) 

C DETERMINE RELATIVE DIFFERENCE 

IF (DIFF .EQ. 0.O) THEN 
RELDIFF = 0.0 

ELSE 

RELDIFF = DIFF/MAX(ABS( VNEXT) , ABSC VTHIS) ) 

END IF 

IF (RELDIFF .LE. EPS) THEN 
MTXKEYC I) = -MTXKEYC I) 

KOUT = KOUT - 1 
END IF 


100 

CONTINUE 



RETURN 


1000 

FORMATC / f f DELETE DUPLICATE VALUES? 

(Y OR NO) : • ,«) 

2000 

FORMATC Al) 


3000 

FORMATC 6X, f ENTER VALUE FOR EPSILON. 

B , R FOR DEFAULT 


* 1PE8.1,’ : 


END 
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SUBROUTINE DELTHRESH( SORTLIST, MTXKEY, ROWKEY, KROW1 , KR0W2 , KOUT, *) 

C 

C DELETE ROWS WHERE THE MERGING TERM IS BELOW A GIVEN THRESHOLD 

C VALUE. 

C 

INTEGER MTXKEYC 200) , ROWKEY(20O) 

REAL S0RTLIST(2, 100) 

CHARACTER* 1 YESNO 
C 

C ASK IF THRESHOLD CULLING IS DESIRED 

TYPE 1000 

ACCEPT 2000 , YESNO 
IF (YESNO .NE. * Y’ ) RETURN 

C SET DEFAULT THRESHOLD 8 REQUEST INPUT VALUE TO REPLACE IT 
THRESH = 1.0 
TYPE 3000, THRESH 
ACCEPT *, THRESH 

C CULL OUT BELOW- THRESHOLD ROWS 

DO 10O 1=1, KR0W1+KR0W2 

IF (MTXKEY(I) .LT. 0) GOTO 100 

IF <SORTLIST(MTXKEY( I) ,ROWKEY( I)) .GE. THRESH) GOTO 150 
MTXKEY(I) = -MTXKEY(I) 

KOUT = KOUT - 1 
100 CONTINUE 

C CHECK TO SEE IF ANY ROWS REMAIN 
IF (KOUT .EQ. 0) THEN 
TYPE *, ’NO MERGING:’ 

TYPE *, ’ALL TERMS ARE LESS THAN THRESHOLD VALUE: ’ , THRESH 
CLOSE( UNIT= 1) 

CLOSE(UNIT=2) 

RETURN 1 
END IF 

C ALL REMAINING TERMS ARE GREATER THAN OR EQUAL TO THRESHOLD 
150 RETURN 

C 

1000 FORMAT( / , ’ DELETE VALUES BELOW A THRESHOLD? (YORN)J ’,*) 

2000 FORMAT( A l ) 

3000 FORMAT* 6X, ’ENTER VALUE FOR THRESHHOLD. ", " FOR DEFAULT ’ , 

* 1PE8.1,’ : ’,«) 

END 
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SUBROUT I NE DELCE I L I NG( SORTL I ST , MTXKEY, ROWKEY, KROW 1 , KR0W2 , KOUT , * ) 


100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

29O0 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3800 

3900 

4000 

4100 


C 

C DELETE ROWS WHERE THE MERGING TERM IS ABOVE A GIVEN CEILING 

C VALUE. 

C 

INTEGER MTXKEY! 200) , ROWKEY( 200): 

REAL SORTL I ST! 2, 100) 

CHARACTER* 1 YESNO 
C 

C ASK IF CEILING CULLING IS DESIRED 

TYPE 1000 

ACCEPT 2000, YESNO 
IF (YESNO .NE. ’ Y’ ) RETURN 

C SET DEFAULT CEILING Q REQUEST INPUT VALUE. TO REPLACE IT 

CEIL = 1.0E+35 
TYPE 3000, CEIL 
ACCEPT *, CEIL 

C CULL OUT ABOVE-CEILING ROWS 

DO 100 1=1, KROW 1 +KR0W2 

IF (MTXKEY(I) .LT. 0) GOTO 100 

IF ( SORTL 1ST ( MTXKEY( I) ,ROWKEY( I)) .GT. CEIL) THEN 
MTXKEY! I) = -MTXKEY(I) 

KOUT = KOUT - 1 
END IF 

100 CONTINUE 

C CHECK TO SEE IF ANY ROWS REMAIN ; 

IF (KOUT .EQ. 0) THEN 

TYPE *, ’NO MERGING: » 

TYPE *, ’ALL TERMS ARE GREATER THAN CEILING VALUE: ' 
CLOSE( UNIT= 1) 

CL0SE(UNIT=2) 

RETURN 1 
END IF 

C ALL REMAINING TERMS ARE LESS THAN OR EQUAL TO CEILING 

150 RETURN 

C 

1000 FORMAT!/,’ DELETE VALUES ABOVE A CEILING? (Y OR N) : ’,«) 

2000 FORMAT! Al) 

3000 FORMAT! 6X, ’ ENTER VALUE FOR CEILING. " , " FOR DEFAULT ’ , 

* 1PE8.1,’ : ’,«) 

END 


CEIL 
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SUBROUTINE GETASCMTX! KROWl , KR0W2 , MROW, NCOL , * ) 
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C 

C OPEN ASSOCIATED MATRIX PAIR (LOGICAL UNIT NOS. 182) AND 
C CHECK FOR VALID DIMENSIONS. 

C 

CHARACTER FILNAM1*35, FILNAM2*35, MATRIX1*8, MATRIX2*8 
C 

C GET NAMES OF ASSOCIATED MATRIX PAIR 
C FIRST MATRIX: 

1 TYPE 1000, ’ 1> ’ , ’ 1’ 

ACCEPT 2000, FILNAM1 

* OPEN(UNIT= 1 , NAME=FILNAM1 , TYPE= ’OLD’ , READONLY, 

* CARR I AGECONTROL= ’LIST’ , ERR= 1 ) 

READ( 1,3000) MATRIX1 , MR0W1 , NC0L1 

C SECOND MATRIX: 

2 TYPE 1000, * 2> ’, * 2 ’ 

ACCEPT 2000, FILNAM2 

OPEN( UNIT=2, NAME=F ILNAM2, TYPE* ’ OLD’ .READONLY, 

* CARR I AGECONTROL= ’ LIST’ ,ERR=2) 

READ( 2 , 3000) MATR I X2 , MR0W2 , NC0L2 

C CHECK FOR VALID DIMENSIONS 

TYPE 4000, ’1> MATRIX1 , MR0W1 , NC0L1 

TYPE 4000, ’ 2> MATR I X2 , MR0V2 , NC0L2 

IF ( ( NCOL 1 . NE . KROW 1 ) .OR. ( NC0L2 . NE . KR0W2 ) ) THEN 

TYPE 5000, ’NO MERGING DONE FOR THIS ASSOCIATED PAIR* 
CLOSE( UNIT= 1 ) 

CLOSE( UN IT=2) 

RETURN 1 

ELSE IF (MR0W1 .NE. MR0W2) THEN 

TYPE 6000, ’NO MERGING DONE FOR THIS ASSOCIATED PAIR* 
CLOSEC UNIT=1) 

CLOSE! UNIT=2) 

RETURN 1 

ELSE 

C MATRICES ARE OK 

MROW = MR0W1 
NCOL = NC0L1 + NC0L2 
RETURN 
END IF 
C 

1000 FORMAT! 1H , A3 , ’ ENTER NAME OF ASSOCIATED MATRIX FILE ’ ,A1,’: ’,«) 
2000 FORMAT! A35 ) 

3000 FORMAT! A8, 2 18) 

4000 FORMAT! 6X, A3 , A8 , ’ IS A ’,18, ’ ROW BY ’,18,’ COLUMN MATRIX’) 

5000 FORMAT! 6X, * # OF COLUMNS DOES NOT MATCH # OF ROWS OF KEY MATRIX’, 

* / , 6X, A40) 

6000 FORMAT! 6X, ’MATRICES HAVE DIFFERENT # OF ROWS’ ,/,6X, A40) 

END 
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100 


206 

C 

300 

C 

400 

C 

300 

600 

700 

C 

800 

C 

900 

C 

1000 

1100 

1200 

1300 

1400 

1 

1500 

C 

1600 

1700 

2 

1800 

1900 

2000 


2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

C 

3200 

C 

3300 

2000 

3400 

3000 

3500 

6000 

3600 

3700 

7000 


SUBROUTINE GETKEYMTX! MATRIX1 , MATRIX2, KROV1 , KR0W2 , KCOL1 ,KC0L2 , *> 

GET KEY MATRIX FILENAMES, OPEN THEM (FORTRAN LOGICAL UNIT 
NOS. 182), READ HEADER RECORDS, AND CHECK FOR COMPATIBLE SIZES. 

CHARACTER*35 FILNAM1 , FILNAM2 
CHARACTERS MATRIX1 , MATRIX2 

FIRST MATRIX: 

TYPE 2000 , ’ 1 * , ’ 1 ’ 

ACCEPT 6000, FILNAM1 

OPEN( UNIT= 1 , NAME= F ILNAM1 , TYPE 3 ’OLD’ .READONLY, 

* CARRIAGECONTROL= ’LIST’ ,ERR= 1) 

READC 1,3000) MATRIX1 , KROWl , KC0L1 

SECOND MATRIX: 

TYPE 2000, ’2’ , ’2’ 

ACCEPT 6000, FILNAM2 

OPEN(UNIT=2,NAME=FILNAM2,TYPE= ’OLD’ .READONLY, 

* CARR I AGECONTROL= ’LIST’ , ERR=2) 

READ(2,3000) MATR I X2 , KROW2 , KC0L2 

CHECK FOR SAME NUMBER OF COLUMNS 
IF (KCOL1 .HE. KC0L2) THEN 
CLOSE! UN IT= 1) 

CLOSE( UNIT=2) 

TYPE 7000, ’ 1> ’, MATRIX1 , KROWl , KCOL1 
TYPE 7000, ’ 2> ’, MATRIX2 , KR0W2 , KC0L2 
TYPE *, ’NO MERGING DONE:’ 

TYPE *, ’MATRICES HAVE DIFFERENT NO. OF COLUMNS’ 

RETURN 1 
END IF 

RETURN 

FORMAT! IX, A 1, ’> ENTER FILENAME OF KEY MATRIX ’ ,A1,’: ’,*) 

FORMAT! A8 , 2 18) 

FORMAT! A35 ) 

FORMAT! 6X, A3 , A8, ’ IS A ’,18,’ ROW BY ’,18,’ MATRIX’) 

END 
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100 

200 

300 

C 

SUBROUTINE GETSORTCOLC MATRIX1 , MATRIX2, KROW1 , KR0W2 , 

* KCOL I , KC0L2 , SORTL 1ST) 


400 

C 

DETERMINE KEY COLUMNS FOR SORTED MERGE, AND READ THEM 


500 

C 

INTO SORTING ARRAY, SORTL I ST. 


600 

700 

800 

C 

REAL S0RTLISTI2, 100) 
CHARACTERS MATRIX1, MATRIX2 


900 

G 


: 

1000 

C 

FIRST MATRIX: 


1100 

1 

TYPE 7000, ’ 1> MATRIX1 , KR0W1 , KC0L1 


1200 


TYPE 8000 


1300 


ACCEPT *, MRGC0L1 

= 

1400 


IF ((MRGC0L1.LT. 1) .OR. ( MRGCOL 1 . GT. KCOL 1 ) ) GOTO 1 

= 

1500 

C 

SECOND MATRIX: 


1600 

2 

TYPE 7000, * 2> ’, MATRIX2 , KR0V2 , KC0L2 

- 

1700 


TYPE 8000 


1800 


ACCEPT * , MRGC0L2 

- 

1900 


IF ( (MRGC0L2.LT. 1) . OR. ( MRGC0L2 . GT. KC0L2) ) GOTO 2 


2000 

C 

READ MERGE COLUMNS INTO SORTING . ARRAY 

- 

2100 


DO 100 J= 1 , MRGC0L1 

- 

2200 


DO 100 1=1, KROW 1 


2300 


READ( 1 , 9000) SORTL IST( 1,1) 

E. 

2400 

100 

CONTINUE 

i 

2500 


DO 200 J= 1 , MRGC0L2 


2600 


DO 200 1=1, KR0W2 

| 

2700 


READ( 2 , 9000) S0RTLIST(2, I) 


2800 

200 

CONTINUE 


2900 


RETURN 

- 

3000 

C 



3100 

7000 

FORMAT! 6X, A3, A8, * IS A ’,18,’ ROW BY ’ , 18, ’ COLUMN MATRIX’) 


3200 

8000 

FORMAT! 9X, ’ENTER COLUMN NO. TO USE FOR SORTED MERGE: ’ ,*) 


3300 

3400 

9000 

FORMAT( E16 . 8) 
END 
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106 


SUBROUTINE KEYL ISTC SORTLIST, KROW1 , KROV2 * MTXKEY, ROWKEY) 


200 

C 



300 

C 

GENERATE KEY LISTS FOR SORTED MERGING 


400 

C 



500 


REAL S0RTLIST(2, 100) 


660 


INTEGER MTXKEY( 100) , ROWKEYt 100) 


700 

C 



800 


11 = 1 


900 


12 = 1 


1000 


DO 100 1=1, KR0W1+KR0W2 


1100 


IF (11 • GT* KROW1) THEN 


1200 

G 

END OF COLUMN FOR MATRIX 1 REACHED; SET KEY FOR 

MATRIX 2 

1300 


MTXKEY(I) = 2 


1400 


ROVKEY(I) = 12 


1500 


12 = 12 + 1 


1600 


ELSE IF (12 .GT. KR0W2) THEN 


1700 

C 

END OF COLUMN FOR MATRIX 2 REACHED; SET KEY FOR 

MATRIX 1 

1800 


MTXKEY(I) = 1 


1900 


ROWKEY(I) = 11 


2000 


11 * 11 + 1 


2100 


ELSE IF (SORTLIST( 1,11) .LE. SORTLIST( 2, 12) ) THEN 


2200 

C 

SET KEY FOR MATRIX 1 


2300 


MTXKEY( I) = 1 


2400 


ROWKEY(I) = 11 


2500 


11 = 11 + 1 


2600 


ELSE 


2700 

c 

SET KEY FOR MATRIX 2 


2800 


MTXKEY(I) = 2 


2900 


ROWKEY(I) = 12 


3000 


12 » 12+1 


3100 


END IF 


3200 

100 

CONTINUE 


3300 


RETURN 


3400 


END 
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SUBROUT I RE MERGEASC ( MTXKEY, KOUT, MROW, NCOL) 

200 

C 


300 

C 

MERGE ASSOCIATED MATRIX PAIR INTO OUTPUT FILE 8 CLOSE FILES 

400 

C 


500 


INTEGER MTXKEY( 200) 

600 


CHARACTER FILNAM3*35, MATRIX3*6 

700 

C 


800 

C 

GET NAME OF RE W ASSOCIATED OUTPUT (MERGED)- MATRIX 

900 


TYPE 1000 

1000 


ACCEPT 2000 , FILRAM3 

1100 


TYPE 3000 

1200 


ACCEPT 4000, MATRIX3 

1300 


OPEH< UR IT= 3 , RAME= FILHAM3 , TYPE= ’ NEW’ , CARRI AGECORTROL= ’ L ISTT ) 

1400 

C 

MERGE MATRIX PAIR INTO OUTPUT MATRIX 

1500 


WRITE! 3, 5000) MATRIX3 , MROW, KOUT 

1600 


DO 100 J= 1 , RCOL 

1700 


DO 100 1=1* MROW 

1800 


READ! ABS< MTXKEYC J) ) * 6000 ) TERM 

1900 


IF (MTXKEY(J) .GT. 0) WRITE( 3,7000) TERM 

2000 

100 

CONTINUE 

2100 

C 

FINISHED: CLOSE FILES 8 RETURN 

2200 


CLOSE! UR I T= 1 ) 

2300 


CLOSE! UNIT=2) 

2400 


CLOSE! UNI T= 3) 

2500 


RETURN 

2600 

C 


2700 

1000 

FORMAT!/, f ENTER NAME OF NEW (MERGED) MATRIX FILE: ’,*) 

2800 

2000 

FORMAT! A35) 

2900 

3000 

FORMAT! 6X, ’ENTER MATRIX NAME: ’ ,0) 

3000 

4000 

FORMAT! A8) 

3100 

5000 

FORMAT! A8, 2 18) 

3200 

6000 

FORMAT! E 16. 8) 

3300 

7000 

FORMAT! 1PE16.8) 

3400 


END 
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SUBROUTINE MERGEKEY< MTXKEY, KOUT) 

MERGE KEY MATRICES INTO OUTPUT MATRIX 8 CLOSE FILES 


iee 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3600 


C 
C 
C 

INTEGER MTXKEYC200) 

CHARACTER'S MATRIX1 , MATRIX2, MATRIX3, FILNAM3*35 
C 

C GET NAME OF OUTPUT KEY MATRIX 8 OPEN NEW FILE 
TYPE 1O00 

ACCEPT 2000, FILNAM3 
TYPE 3000 

ACCEPT 4000, MATRIX3 

OPEN ( UN I T= 3 , NAME= F I LNAM3 , TYPE= * NEW’ , C ARRI AGECONTROL- * LIST V ) 
C MERGE KEY MATRICES INTO OUTPUT MATRIX 
REWIND 1 
REWIND 2 

READ( 1 , 5000 ) MATRIX! ,XR0W1 , KCOL 1 
READ( 2,5000) MATRIX2 , KR0W2 , KC0L2 
WRITE( 3 , 5000) MATRIX3, KOUT, KCOL 1 
DO 100 J= 1 , KC0L1 
DO 100 1=1, KROW1+KROW2 

READ( ABS( MTXKEYC I)) , 6000 ) TERM 
IF (MTXKEY(I) .GT. 0) WRITE( 3, 7000) TERM 
100 CONTINUE 

C FINISHED: CLOSE FILES 8 RETURN 
CLOSE( UNIT= 1 ) 

CLOSE( UNIT=2) 

CL0SE(UNIT=3) 

RETURN 

C 

1000 FORMATC / , * ENTER NAME OF NEW (MERGED) MATRIX FILE: ',#) 

2000 FORMAT( A35) 

3000 FORMATC 6X, * ENTER MATRIX NAME: * ,$) 

4000 FORMATC A8) 

5000 FORMATC AS, 2 18) 

6000 FORMATC El 6. 8) 

7000 FORMATC 1PE16 . 8) 

END 
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STATCQRR 


The statistical correlation program STATCORR was written for implementation on the GSFC 
Code 750 VAX-1 1/780 computer system, and the analysis was done using the Code 750 VAX ver- 
sion of Level 17.5 NASTRAN. The empirical testing was performed by the GSFC Environmental 
Test & Integration Branch and test mode shapes, frequencies, and damping were determined us- 
ing the standard in-house processing software on a PDP-1 1/35 computer system. 

Three data files are required to input analytical mode shape data to the STATCORR pro- 
gram: 1) the LAMA matrix file, which contains data from the NASTRAN LAMA table ( frequen- 
cy ; mass, and stiffness vectors) and mode symmetry information ; 2) the PHITE matrix file, which 
contains the eigenvectors generated by the NASTRAN run; 3) the grid point list file, obtained 
from the TESET vector. All the data required to build these files is contained in the NASTRAN 
punch file produced by using the DMAP ALTER package discussed earlier. 

Since the NASTRAN analysis and the statistical correlation are performed on the same com- 
puter system, the NASTRAN punch file is immediately available for processing (to generate the 
analytical input files to STATCORR) by the appropriate intermediate processor programs 
MERGE, LAMA, UNPACKDMI, GRDPTLIST, discussed subsequently. These processor programs 
read the NASTRAN punch file and produce input files in the format required for input to 
STATCORR. 

One data file for each experimental mode is required to input experimental mode shape data 
to the STATCORR program. Each experimental data file contains all data (frequency, damping, 
symmetry, gridpoint IDs, and mode shape data) to describe one experimental mode. The in-house 
Modal Survey processing program has the capability of generating various ASCII data files contain- 
ing the necessary grid-point ID, frequency, damping, and mode-shape information, which can be 
further processed to generate files in the appropriate format for input to the STATCORR program. 

Since the experimental mode data is stored on a different computer system, it was necessary 
to provide a method of data transfer between the Code 750 VAX-1 1/780 computer and the Envi- 
ronmental Test & Integration Branch PDP-11/35 computer. Because rapid, convenient file trans- 
fer was desired and the two computer systems had no similar storage media to permit simple vol- 
ume transfer between the machines, a remote terminal emulator program with ASCII file transfer 
capability was adapted for use on the PDP-1 1/35. This permits the PDP computer to be used as 
a remote (dial-up) terminal to the VAX-1 1/780, allowing direct transfer of ASCII text files via 
phone line between the VAX and PDP computers. 

For a user’s guide describing the execution of NASTRAN, STATCORR, and the associated 
preprocessing programs, please refer to NASA Technical Memorandum 86044. Also included is 
a working example using the SPARTAN-1 model, and FORTRAN source code listings of the 
STATCORR program and the preprocessing programs TESETDMI, UNPACKDMI, LAMA, and 
GRDPTLST. 
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* * * * 


PROGRAM STATCORR 
C 

PARAMETER (MAXMODE= 200) 
PARAMETER ( MAXDOF = 1000) 

COMMON /LIMITS/ MXMD.MXDF 
DATA MXMD/MAXMODE/ , MXDF/MAXDOF/ 


REAL FREQAN! MAXMODE) , FREQEX! MAXMODE) 

REAL DISPANI MAXDOF, MAXMODE) , DISPEXI MAXDOF, MAXMODE) 

REAL MASS< MAXMODE) , STIFF! MAXMODE) , DAMP! MAXMODE) 

REAL CORREL! MAXMODE, MAXMODE) , C( MAXMODE, MAXMODE) 

REAL S( MAXMODE, MAXMODE) 

REAL RMSA( MAXMODE, MAXMODE) , RMSE< MAXMODE, MAXMODE) 

REAL AVECT( MAXDOF) , EVECT( MAXDOF) , DIFF( MAXDOF) , DIFABSI MAXDOF) 
INTEGER IDDOFA! MAXDOF, MAXMODE) , I DDOFEI MAXDOF, MAXMODE) 

INTEGER ICOMPA! MAXDOF, MAXMODE) , ICOMPE! MAXDOF, MAXMODE) 
INTEGER NDOFA( MAXMODE) , NDOFE( MAXMODE) 

INTEGER IDDOFK MAXDOF) , ICOMP I (MAXDOF) 

INTEGER NAF IT( MAXMODE) , MEF IT( MAXMODE) 

I NTEGER ASYMM( 3 , MAXMODE) , ESYMM< 3 , MAXMODE) 

CHARACTER* 10 AHEAD ER( 4) , EHEADER( 4) , HEADER! 32) 

CHARACTERS CDATE 
CHARACTER* 1 YESNO 
LOGICAL PRTAN, PRTEX, LSYMM, LPRINT 

COMMON /WHEN/ CDATE 


C 


C 

c 


c 


c 


c 


c 

c 

c 


c 

c 


DATA PRTAN , PRTEX, LSYMM, LPR I NT/4* . FALSE . / 

DATA AIIEADER/2*’ NASTRAN’ ,2*' TEST’/ 

DATA EHEADER/2* * TEST’, 2*’ NASTRAN’/ 

DATA HEADER/ ’ * OF GP’’S’, ’ CORREL’ , 2*’ 2*’ RMS’, 

’ MAX REL’ , ’ GRID’, ’ MODE’, ’FREQUENCY’, 

’ MODE’, ’FREQUENCY’, ’ COMPARED’, * COEFF’ , 
’ CA’, ’ S’, ’(NASTRAN)’, ’ (TEST)’, 

’DIFFERENCE’, ’ POINT’, 12*’ .’/ 


PRINT SECTION 1 HEADER 
CALL DATE( CDATE) 
PRINT 80G0, CDATE 
PRINT 4000 


CHECK FOR SEPARATE PRINT FILE 
TYPE 5000 

ACCEPT 7000, YESNO 

IF (YESNO .EQ. ’Y’) LPRINT = .TRUE. 

IF (LPRINT) PRINT 5500, YESNO 

GET ANALYTICAL FREQ, MASS, STIFFNESS, MODE-SHAPE, 3 DOF LISTS 
CALL GET AN ( FREQAN, MASS, STIFF, AS YMM, DISPAN, NAN, 

* * IDDOFA, ICOMPA.NDOFA, LPRINT) 

GET EXPERIMENTAL FREQ, DAMPING, MODE-SHAPE, 8 DOF LISTS 

CALL GETEXP( FREQEX, DAMP, ESYHM, DISPEX.MEX, IDDOFE, ICOMPE, 

^ NDOFE LPRINT) 

DETERMINE SYMMETRY, MODE-SHAPE DUMP, 8 THRESHOLD OPTIONS 
CALL OPTIONS! LSYMM, PRTAN, PRTEX, RTHRESH, LPRINT) 

PRINT INPUT SUMMARY 

CALL INPSUM(NAN, FREQAN, MASS, STIFF, ASYMM, MEX, FREQEX, DAMP, ESYMM) 
SORT DOF LISTS 3 MODE SHAPE MATRICES (BASED ON DOF SORT) 

CALL DOFSORT! IDDOFA, ICOMPA.NDOFA, DISPAN, NAN) 

CALL DOFSORT! IDDOFE, ICOMPE, NBOFE.B I SPEX, MEX) 

COMPUTE CORRELATION COEFFICIENTS 8 RMS VALUES FOR ALL 

POSSIBLE ANALYTICAL/ EXPERIMENTAL PAIRS 
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G 


G 


100 

c 

c 

c 

c 


c 


c 


c 


200 

C 


c 


c 


c 


300 

C 


DO 100 R= 1 , NAN 
DO 100 M= 1 , MEX 

GET VECTORS TO BE COMPARED 0 THEIR DOF INTERSECTION i SET 
CALL GETVEC! D I SPAN! 1 , N) , NDOFA< N) , AVECT) 

CALL GETVEC( D I SPEX( 1 , M) , NDOFE( M) , EVECT) 

CALL INTERSECT! AVECT, EVECT, I DDOF A( 1 , N) , IDDOFE! 1 , M) , 

* ICOMPA! 1 , N) , ICOMPE! 1 , M) , 

* NDOFA( N> , NDOFE! M) , I DDOF I , ICOMP I , NDOF I , 

* ASYMM( 1 , N) , ESYMM< 1 , M) , LSYMM) 

COMPUTE CORR COEFFS 0 RMS VALUES FOR THIS PAIR 

CALL CORRMS ( AVECT , EVECT , NDOF I , CORREL! N , M) , 

* C(N,M) ,S(N,M) , RMSA( N,M) ,RMSE( N, M) ) 


CONTINUE 

PRINT CORRELATION COEFFICIENT TABLE 
CALL CORRTBL! CORREL , NAN , MEX) 

DETERMINE BEST MATCH FOR EACH EXPERIMENTAL 0 ANALYTICAL MODE 
CALL MATCH! CORREL, NAN, MEX, NAF IT, MEFIT) 

FOR EACH MATCHED PAIR, GET X/RMS DIFFERENCES 0 PRINT SUMMARY 
PRINT ANALYTICAL MODE SHAPE HEADER 
PRINT 8000, CDATE 
PRINT 1000, AHEADER, HEADER 
DO 200 N= 1 , NAN 
M = NAF IT! N) 

GET VECTORS TO BE COMPARED 0 THEIR DOF INTERSECTION SET 
CALL GETVEC! DISPAN! 1 , N) , NDOF A! N) , AVECT) 

CALL GETVEC! D I SPEX! 1 , M) , NDOFE! M) , EVECT) 

CALL INTERSECT! AVECT, EVECT, I DDOF A! 1 , N) , IDDOFE! 1 , M) , 

* ICOMPA! I . N) , ICOMPE! 1 , M) , 

* HDOFA! N) , NDOFE! M) , I DDOF I , ICOMP I , NDOF I , 

* ASYMM! 1 , N) , ESYMM! 1 , II) , LSYMM) 

CALCULATE INDIVIDUAL DIFFERENCES 

CALL RMSDIFF! AVECT, EVECT, NDOFI, RMSA! N,M) .RMSE! N,M) , 

* CORREL! N, M) , S! N, M) , DIFF , DIFABS) 

PRINT SUMMARY FOR THIS PAIR, WITH MAX DIFF 0 > THRESHOLD 
C ALL PRINT! N,FREOAN!N) , M, FREQEX! M) .NDOFI .CORREL! N,M) , 

* C! H,M) ,S!N,M) , RMSA! N , M) , RMSE! N , M) , 

« DIFF , DIFABS , IDDOFI , ICOMP I , RTHRESH) 


CONTINUE 

PRINT EXPERIMENTAL MODE SHAPE HEADER 
PRINT 8000, CDATE 
PRINT 3000, EHEADER, HEADER 
DO 300 M= 1 , MEX 
N = MEFITIM) 

GET VECTORS TO BE COMPARED 8 THEIR DOF INTERSECTION SET 
CALL GETVEC! D I SPEX! 1 , M) , NDOFE! M) , EVECT) 

CALL GETVEC! DISPAN! 1 , N) , NDOFA! N) , AVECT) 

CALL INTERSECT! AVECT, EVECT, IDDOFA! 1 , N) , IDDOFE! 1 , M) , 

* ICOMPA! 1 , N) , ICOMPE! 1 , M) , 

* NDOFA! N ) , NDOFE! M) , I DDOF I , I COMP I , NDOF I , 

* ASYMM! 1 , N) , ESYMM! 1 , M) , LSYMM) 

CALCULATE COMBINED RMS AND INDIVIDUAL DIFFERENCES 

CALL RMSDIFF! AVECT, EVECT, NDOF I, RMS A! N,M) , RMSE! N , M) , 

* CORREL!N,M) ,S!N,M) .DIFF, DIFABS) 

PRINT SUMMARY FOR THIS PAIR, WITH MAX DIFF 0 > THRESHOLD 
CALL PRINT! M.FREQEXIM) , N, FREOAN! N) , NDOFI .CORREL! N, M) , 

* C! N, M) ,S! N, M) , RMSA! N, M) .RMSE! N, M) , 

* DIFF, DIFABS, IDDOFI, I COMP I, RTHRESH) 


CONTINUE 

PRINT ANALYTICAL 0/OR EXPERIMENTAL MODE-SHAPE VECTORS 
CALL PRINTPHI! NAN, NDOFA, IDDOFA, ICOMPA.DISPAN.PRTAN, 
* MEX, NDOFE, IDDOFE, ICOMPE, D I SPEX, PRTEX) 
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STOP ’FORTRAN STOP — PROCESSING COMPLETED’ 

C 

1000 FORMAT( 1 H0 , * 4 . ANALYTICAL MODE SHAPES AND THEIR BEST’ * 

* ’ EXPERIMENTAL MATCHES: ’,/’ 0’ , 12( A10, IX) ,2! /12< IX, A10) ) ) 
3000 FORMAT! 1H0, ’5. EXPERIMENTAL MODE SHAPES AND THEIR BEST’ 

* ’ ANALYTICAL MATCHES :’,/’ 0 ’,12! A10 , IX) , 2! / 12( IX, A10) ) ) 

4000 FORMAT( 1H0, ’ 1 . INTERACTIVE DIALOG: • ) 

3000 FORMAT( /, ’ IS A SEPARATE OUTPUT LISTING FILE TO BE PRINTED?’ , 

* ’ ( Y OR N) : ’ ,6) 

5500 FORMAT!/,’ IS A SEPARATE OUTPUT LISTING FILE TO BE PRINTED?’, 

* ’(Y OR N) : ’ ,A1) 

7000 FORMAT! Al) 

8000 FORMAT! 1 HI, ’ NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL’, 

* ’ CORRELATION’, /IX, A9) 

END 
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SUBROUTINE CORRMS ( AVECT, EVECT , NDOF I , CORREL , CA, S , RMS A, RMSE) 

COMPUTE CORRELATION COEFFICIENTS 8 RMS VALUES 

REAL AVECT(*) ,EVECT(*) 

INITIALIZE RMS’S 8 CORE COEFFS 
RMSA = 0*0 
RMSE = 0.0 
CORREL =0.0 
CA = 0.0 
S = 0.0 

IF (NDOF I .EQ. 0) RETURN 

COMPUTE VARIANCES OF ANALYTICAL 8 EXPERIMENTAL VECTORS 
NOTE: THE FINAL VALUES COMPUTED HERE FOR THE "VARIANCES B 
ARE NOT TRUE VARIANCES. TO OBTAIN VARIANCES, THE VALUES 
VARA, VARE, 8 VARAE SHOULD BE EACH MULTIPLIED BY 1/NDOFI . 
THESE VALUES ARE NOT REPORTED, BUT ARE USED ONLY TO 
CALCULATE CORREL, CA, S, RMSA, 8 RMSE. THE 1/NDOFI FACTOR 
IS NOT INCLUDED, SINCE IT APPEARS IN BOTH THE NUMERATOR AND 
DENOMINATOR OF CORREL 8 CA; THE EXTRA ARITHMETIC OPERATIONS 
WOULD SIMPLY INTRODUCE UNNECESSARY ERROR IN THE RESULTS. 

THE 1/NDOFI FACTOR IS INCLUDED IN THE CALCULATIONS OF 
S, RMSA, 8 RMSE. 

VARA = 0.0 
VARE = 0.0 
VARAE = 0.0 
DO 100 N= 1 , NDOF I 

VARA = VARA + AVECTC N) *AVECT( N) 

VARE = VARE + EVECT( N) *EVECT( N) 

VARAE = VARAE + AVECTC N) *EVECT< N) 

CONTINUE 

COMPUTE CORRELATION COEFFICIENT FOR THIS PAIR OF VECTORS 
IF (VARAE .NE. 0.0) CORREL = VARAE/SQRT( VARA* VARE) 
COMPUTE CORREL COEFF (REFERENCED TO ANALYTICAL VECTOR) 

IF (VARA .NE. 0.0) CA = VARAE/VARA 
COMPUTE ROOT OF MEAN SQUARE DIFFERENCE 

S = SGRT( ABS( VARE- C A* C A* VARA) /NDOF I ) 

COMPUTE RMS VALUES 

NOTE: FOR THIS APPLICATION, RMS = STD. DEV. = SQRT( VARIANCE) 

RMSA = SQRT( VARA/NDOFI) 

RMSE = SQRT( VARE/NDOFI) 

RETURN 

END 
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SUBROUTINE CORRTBU CORREL, NAN, MEX) 

C 

C PRINT AN NAN X MEX TABLE OF THE ANALYTICAL VS. EXPERIMENTAL 
C CORRELATION COEFFICIENTS. 

C 

COMMON /LIMITS/ MAXMODE , MAXDOF 

COMMON /WHEN/ CDATE 

REAL CORREL ( MAXMODE, MAXMODE) 

CHARACTER CDATE*9 
C 

DO 200 MFIRST= 1 , MEX t 16 
DO 200 NF IRST= 1 , NAN , 50 

ML AST = MIN( MEX. MFIRST+15 ) 

NLAST = MIN( NAN, NFIRST+49 ) 

C PRINT HEADER FOR ONE PAGE OF THE TABLE 

PRINT 1000, CDATE 
PRINT 2000 

PRINT 3000, ( M, M=MFIRST, MLAST) 

DO 100 N=NF IRST, NLAST 

PRINT 4000, N, (CORREL(N,M) ,M=MFIRST,MLAST) 

100 END DO 

200 END DO 

RETURN 

1000 FORMAT( 1H1 , ’ NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 

* ’CORRELATION’, /1X,A9, 

* / 1 HO , ’ 3 . CORRELATION COEFFICIENTS FOR ANALYTICAL VS. 

* ’EXPERIMENTAL COMPARISONS: ’ ) 

2000 FORMAT< ’©ANALYTICAL’ , 24X, ’EXPERIMENTAL MODES’) 

3000 F0RMAT(3X, ’MODES’ ,3X, 15, 1517) 

4000 FORMAT( IX, 16, 4X, 16( 1X.F6.3) ) 

END 
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anna 


SUBROOT I RE DOFSORT( I DDOF , I COMP , RDOF , D ISP , RHODES) 

PUT DOF LIST FOR EACH MODE IRTO ASCERDIRG SORT, 

ARD SLAVE-SORT THE CORRESPORDIRG MODE-SHAPE VECTORS. 


C 


100 

200 


COMMOR /LIMITS/ MAXMODE , MAXDOF 
REAL D ISP (MAXDOF, MAXMODE) 

I RTEGER I DDOF ( MAXDOF , MAXMODE) , I COMP ( MAXDOF , 


MAXKODE) , RDOF ( MAXMODE) 


DO 


200 

IF 

DO 

DO 


J= 1 , RMODES 
(RDOF(J) .EQ. 


1) GOTO 200 


100 1= 1,HD0F( J)-l 
100 K= 1+ 1 , RDOF( J) 

IF (IDDOF(K.J) .LT. IDDOF(I,J)) THEM 
CALL SVAPC I DDOF ( I , J ) , IDDOF(K.J) ) 

CALL SWAP( ICOMP(I,J), ICOMP(K.J) ) 

CALL SWAP ( DISP(I,J), DISP(K,J) ) 

ELSE IF ((IDDOF(K.J) .Eft. IDDOF(I,J)) .ARD. 
(ICOMP(K.J) .LT. ICOMP(I.J)) 

SWAP( IDDOF(I,J), IDDOF(K.J) ) 
SWAP( ICOMP(I.J), ICOMP(K.J) ) 
SWAP< DISP(I,J), DISP(K.J) ) 


CALL 

CALL 

CALL 

IF 


END 
CONTINUE 
CONTINUE 
RETURN 
END 


THEN 
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c 

c 

c 

c 


SUBROUT I NE GETAN ( FREQS , MASS .STIFF, AS YMM; SHAPES , NMODES , 

* IDDOF , ICOMPA, NDOFS . LPRINT) 

GET FREQUENCY, MASS, STIFFNESS 8 SYMMETRY VECTORS; 

MODE-SHAPE MATRIX; AND GRID POINT LIST. 

COMMON /LIMITS/ MAXMODE , MAXDOF 


C 

C 

ie 


20 


30 


c 


40 


REAL FREQS ( MAXMODE) , MASS< MAXMODE) , STIFFC MAXMODE) 
REAL SHAPES( MAXDOF, MAXMODE) 

INTEGER ASYMM( 3, MAXMODE) 

INTEGER IDDOF ( MAXDOF , MAXMODE) , ICOMPA( MAXDOF , MAXMODE) 
INTEGER NDOFS( MAXMODE) 

CHARACTER*35 FRQFIL,SHPFIL,DOFFIL 
CHARACTER'S FRQMTX, SHPMTX 
LOGICAL LPRINT 


GET FILE NAMES 
TYPE 1000 

ACCEPT 2000, FRQFIL 

OPEN(UNIT= 1,NAME=FRQFIL,TYPE= ’OLD* , READONLY, ERR= 10) 

IF (LPRINT) PRINT 1500, FRQFIL 
TYPE 3000 

ACCEPT 2000, SHPFIL 

OPEN( UNIT=2, NAME=SHPFIL, TYPE= ’OLD* , READONLY, ERR= 20) 

IF (LPRINT) PRINT 3500, SHPFIL 
TYPE 4000 

ACCEPT 2000, DOFF I L 

OPEN( UN I T= 3 , NAME= DOFF I L , TYPE= ' OLD ’ , READONLY, ERR= 30 ) 

IF (LPRINT) PRINT 4500, DOFFIL 
READ FREQUENCY FILE 

READ( 1 , 7000) FRQMTX, NFREQS , NCOLS 
IF (NCOLS .NE. 6) THEN 

TYPE 5000, FRQMTX, NCOLS 

IF (LPRINT) PRINT 5000, FRQMTX, NCOLS 

CLOSE( UNIT 3 1 ) 

CL0SE(UNIT=2) 

CLOSE! UN I T= 3) 

STOP 
END IF 

IF (NFREQS .GT. MAXMODE) PRINT 6000, NFREQS, MAXMODE 
NMODES = M I N( NFREQS, MAXMODE) 

DO 100 1=1,6 
DO 100 J=l, NFREQS 

IF (J .GT. MAXMODE) THEN 
READ( 1,8000) DUMMY 

ELSE IF (I .EQ. 1) THEN t GET FREQUENCY LIST 
READ( 1,8000) FREQS(J) 

ELSE IF (I .EQ. 2) THEN ! GET MASS LIST 
READ( 1,8000) MASS(J) 

ELSE IF (I .EQ. 3) THEN ! GET STIFFNESS LIST 
READ( 1,8000) STIFF(J) 

ELSE IF (I .EQ. 4) THEN I GET PLANE 1 SYMMETRY LIST 
READ( 1,8000) SYMM 
ASYMM(l.J) = NINT(SYMM) 

else r1Id( ( !,80 E ?0) 5 ^ n r GET PLANE 2 SYMMETRY LIST 

ASYMM( 2 , J) = NINT(SYMM) 

ELSE IF (I .EQ. 6) THEN ! GET PLANE 3 SYMMETRY LIST 
READ( 1,8000) SYMM 
ASYMM( 3 , J) = NINT(SYMM) 
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END IF 

100 CONTINUE 

C READ MODE-SHAPE FILE 

READ( 2,7000) SHPMTX, NDISPS , NSHAPES 
IF (NSHAPES . NE. N FREQS) THEN 

NMODES = MIN( NFREQS, NSHAPES , NMODES ) 

TYPE 9000 , NFREQS , NSHAPES , NMODES 

IF ( LPRINT) PRINT 9000, NFREQS , NSHAPES , NMODES 

END IF 

IF (NDISPS . GT. MAXDOF) PRINT 10000, NDISPS, MAXDOF 
DO 200 J=l, NMODES 
DO 200 1=1, NDISPS 

IF (I . LE. MAXDOF) READ( 2,8000) SHAPES! I, J> 

IF (I . GT. MAXDOF) READ(2,8000) DUMMY 
200 CONTINUE 

C READ DOF ID FILE 

NDOFS(l) = 0 
DO 300 1= 1, 1000000 

IF (I . LE. MAXDOF) 

* READ( 3, 1 1 000 , END= 350) IDDOF( 1,1), ICOMPA( 1,1) 

IF (I .GT. MAXDOF) READ! 3, 1 1000, END=350) IDUMMY, IDUMMC 
NDOFS(l) = NDOFS(l) + 1 
300 CONTINUE 

350 IF (NDOFS(l) .NE. NDISPS) THEN 

TYPE 12000, NDISPS, NDOFS( 1) 

IF (LPRINT) PRINT 12000, NDISPS, NDOFS( 1 ) 

CLOSE! UNIT= 1 ) 

CLOSE! UNIT=2) 

CLOSE! UNIT 2 3) 

STOP 
END IF 

NDOFS(l) = MIN(NDOFS( 1) , MAXDOF) 

IF (NMODES .EQ. 1) RETURN 
DO 500 J= 2, NMODES 

NDOFS(J) = NDOFS(l) 

DO 400 1=1, NDOFS! J) 

I DDOF ! I , J ) * IDDOF! 1,1) 

ICOMPA! I , J) = ICOMPA(I,l) 

400 CONTINUE 

500 CONTINUE 

CLOSE! UN I T= 1) 

CLOSE! UN I T= 2) 

CLOSE! UN I T= 3) 

RETURN 

C 

1000 FORMAT!/, • ENTER ANALYTICAL LAMA MATRIX FILENAME: 1 ,0) 

1500 FORMAT! / , * ENTER ANALYTICAL LAMA MATRIX FILENAME: ’,A35) 

2000 FORMAT! A35) 

3000 FORMAT! 1 ENTER ANALYTICAL MODE-SHAPE MATRIX FILENAME: * ,0) 

3500 FORMAT! 1 ENTER ANALYTICAL MODE-SHAPE MATRIX FILENAME: *,A35) 

4000 FORMAT! ’ ENTER ANALYTICAL GRID POINT LIST FILENAME: \S) 

4500 FORMAT! 1 ENTER ANALYTICAL GRID POINT LIST FILENAME: \A35) 

5000 FORMAT! 1 FREQUENCY MATRIX ’,A8,’ HAS ’,16,’ COLUMNS. 1 

* 1 SHOULD BE 6 COLUMNS . 1 ) 

6000 FORMAT! 1 **** WARNING: FREQUENCY VECTOR HAS 1 , 13, 1 ENTRIES. 1 , 

* /15X, ’ONLY THE FIRST 1 , 15, 1 WILL BE USED. 1 ) 

7000 FORMAT! A8, 2 18) 

8000 FORMAT! El 6. 8) 

9000 FORMAT! 1 **** WARNING: UNEQUAL NUMBER OF FREQUENCIES AND 1 

* 1 MODE SHAPES . 1 , 

* /15X, ’NUMBER OF FREQUENCIES: 1 ,^, 


49 



* / 15X, * NUMBER OF MODE SHAPES: ’, 16, 

* /15X, ’ONLY THE FIRST’ ,13,’ WILL BE USED.’) 

10000 FORMAT! ’ **** WARNING: MODE-SHAPE VECTORS HAVE ’,13,’ ENTRIES. ’ , 

* /15X, ’ONLY THE FIRST’, 16,’ WILL BE USED.’) 

1 1000 FORMAT! 18, IX, 1 1) 

12000 FORMAT! ’ **** ERROR: UNEQUAL NUMBER OF MODE-SHAPE POINTS AND ’ , 

* ’D.O.F. ID’ 'S. ’,16,’ 8 ’,16) 

END 
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c 

c 

c 

G 


C 

C 


G 

ie 


c 


c 

100 

150 

200 

250 


C 

1000 

1500 

2O00 

5000 

6000 

END 


SUBROUTINE GETEXP( FREQEX, DAMP y ESYMM, DISPEX, MEX, 

* IDDOFE, I COMP , NDOFE, LPRINTJ 

GET EXPERIMENTAL FREQUENCY 8 DAMPING LISTS, MODE-SHAPE MATRIX, 
AND D.O.F. LIST. 

COMMON /LIMITS/ MAXMODE , MAXDOF 

REAL FREQEX( MAXMODE) , DAMP! MAXMODE) , DISPEX! MAXDOF, MAXMODE) 

REAL DUMMY ( 5) 

INTEGER ESYMM( 3 , MAXMODE) 

INTEGER IDD0FE( MAXDOF , MAXMODE) , ICOMP< MAXDOF, MAXMODE) 

INTEGER NDOFE( MAXMODE) 

CHARACTER* 35 MSFILE 
LOGICAL LPRINT 

PROCESS MODE-SHAPE FILES 

TYPE *, • ’ f SKIP A LINE 
IF (LPRINT) PRINT *, ’ 1 
MEX = 0 

DO 200 J= 1 , MAXKODE 

GET NEXT EXPERIMENTAL MODE-SHAPE FILE 
TYPE 1000 

ACCEPT 2000, MSFILE 

IF (MSFILE .EQ. * NONE* ) GOTO 250 

OPEN( UN IT= 10 , NAME=MSF ILE, TYPE= ’ OLD’ , READONLY, ERR= 10) 
IF (LPRINT) PRINT 1500, MSFILE 
GET FREQUENCY, DAMPING 8 SYMMETRY FOR THIS MODE-SHAPE 
READ( 10 , *) FREQEX(J) 

READ( 10 , *) DAMP(J) 

READ( 10 , *) ( ESYMM( K, J) , K= 1 , 3) 

GET MODE-SHAPE DISPLACEMENTS 8 D.O.F. ’S 
DO 100 1=1, MAXDOF 

READ( 10, * , END= 150) 

* IDDOFE( I, J) , ICOMP( I, J) , DISPEXC I, J) 
NDOFE(J) = I 
CONTINUE 

PRINT 5000, MAXDOF, MAXDOF 
MEX = J 
CLOSE( UNIT= 10) 

CONTINUE 

TYPE 6000, MAXMODE 
IF (LPRINT) PRINT 6000, MAXMODE 
IF (LPRINT .AND. (MSFILE . EQ. ’NONE’)) PRINT 1500, 

IF (MEX .EQ. 0) THEN 

TYPE **** ERROR: NO EXPERIMENTAL MODE-SHAPE 
IF (LPRINT) PRINT *, 

* • * *** ERROR: NO EXPERIMENTAL MODE-SHAPE 
STOP 
END IF 

RETURN 


MSFILE 
REQUESTS * 
REQUESTS ’ 


* 

* 


* 

* 


FORMAT( ’ ENTER NEXT EXPERIMENTAL MODE-SHAPE FILENAME. ’ , 
’("NONE" IF NO MORE) : ’,$) 

FORMAT( * ENTER NEXT EXPERIMENTAL MODE-SHAPE FILENAME. 

’ ( " NONE " IF NO MORE) : 1 ,A35) 

FORMAT( A35 ) 

FORMAT( * **$* WARNING: MODE-SHAPE VECTOR HAS*,I6,* OR MORE’, 

1 ENTRIES. ’ ,/13X, ’ONLY THE FIRST’ , 16, ' WILL BE USED. ’ ) 
FORMAT! * **** NO MORE EXPERIMENTAL MODE-SHAPE FILES PERMITTED. 
’ MAXIMUM IS • , 14) 
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SUBROUTINE GETVEC( D ISP, NDOF, VECT) 


C 

100 


COPY VECTOR OF LENGTH NDOF FROM PERMANENT VECTOR TO 
WORKING VECTOR. 

REAL DISP( *) , VECT( *) 

IF (NDOF . EQ. 0) RETURN 
DO 100 1=1, NDOF 

VECT(I) = DISP(I) 

CONTINUE 

RETURN 

END 
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SUBROUTINE INPSUM! NAN, FREQAN, MASS, STIFF, ASYMM, 

* MEX, FREQEX, DAMP , ESYMM) 

C 

C PRINT A SUMMARY OF THE ANALYTICAL 8 EXPERIMENTAL FREQUENCIES, 

C MASS, STIFFNESS, DAMPING, AND SYMMETRY VECTORS. 

C 

REAL FREQAN!*) .MASS!*) , STIFF!*) 

REAL FREQEX! *) , DAMP( #) 

INTEGER ASYMM(3,*) ,ESYMM!3,*) 

LOGICAL NOSYMAC 3) , YESYMA! 3) ,NOSYME<3) , YESYME! 3) .FIRST 
CHARACTER CDATE*9 
COMMON /WHEN/ CDATE 

DATA NOSYMA, YESYMA.NOSYME, YES YME/ 12*. FALSE./, FIRST/. TRUE./ 

C 

C PRINT HEADER 

PRINT 8000, CDATE 
PRINT 1000 

C PRINT ANALYTICAL SUMMARY 

C PRINT ANALYTICAL HEADER 

PRINT 2000 

C PRINT SUMMARY DATA 

PRINT 3000, ( I,FREQAN( I) ,MASS( I) .STIFF! I) , 

* ( ASYMM! J , I ) , J= 1 , 3) , 1=1, NAN) 

C CHECK FOR INCONSISTENT ANALYTICAL "UNDEFINED" SYMMETRY 

DO 200 IAXS= 1 , 3 

DO 100 MODEM, NAN 

IF ! ASYMM! I AXS, MODE) .EQ.0) NOSYMA! IAXS) = . TRUE. 
IF ! ASYMM! IAXS, MODE) .NE.0) YESYMA! IAXS) = .TRUE. 
100 CONTINUE 

IF ! NOSYMA! IAXS) .AND. YESYMA! IAXS)) PRINT 6000, IAXS 
200 CONTINUE 

C PRINT EXPERIMENTAL SUMMARY 

C PRINT EXPERIMENTAL HEADER 

PRINT 4000 

C PRINT SUMMARY DATA 

PRINT 5000, ! I.FREQEX! I) .DAMP! I) , 

* ! ESYMM! J , I) , J= 1 , 3) , 1=1, MEX) 

C CHECK FOR INCONSISTENT EXPERIMENTAL "UNDEFINED" SYMMETRY 

DO 400 I AXS= 1 , 3 

DO 300 MODEM, MEX 

IF ! ESYMM! IAXS, MODE) .EQ.0) NOSYME! IAXS) = .TRUE. 
IF ! ESYMM! I AXS, MODE) .NE.0) YES YME! IAXS) =. TRUE. 
300 CONTINUE 

IF ! NOSYME! IAXS) .AND. YESYME! IAXS) ) PRINT 6000, IAXS 
400 CONTINUE 

C CHECK FOR INCONSISTENT ANA VS. EXP "UNDEFINED" SYMMETRY 

DO 500 IAXSM.3 

IF ! ! NOSYMA! IAXS) .AND. YESYME! IAXS)) .OR. 

* ! YESYMA! IAXS) . AND. NOSYME! IAXS) ) ) THEN 

IF ! FIRST) PRINT 7000 

FIRST = .FALSE. 

PRINT 6000, IAXS 
END IF 

500 CONTINUE 

RETURN 
C 

1000 FORMAT! 1H0, '2. SUMMARY OF FREQUENCY, MASS, STIFFNESS,’, 

* ’ DAMPING, 8 SYMMETRY:’) 

2000 FORMAT! // , 5X, ’ ANALYTICAL MODES : ’ 

* /1H0.4X, ’MODE ’,’ FREQUENCY ’,’ MASS 

* ’ STIFFNESS ’ . ’ SYMMETRY 
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* / , 5X, * ’ ,4( IX, ’ ’)) 

3000 FORMAT! <NANX/5X, 13, IX, 3( 1PE16.8) ,3X,3I3) ,/) 

4000 FORMAT! // , 5X, * EXPERIMENTAL MODES : f , 

* / 1 H0 , 4X, 1 MODE FREQUENCY DAMPING SYMMETRY’ , 

^ / 5X ’ ’ 3( 1 X ’ * ) ) 

5000 FORMAT( < MEX> ( /5X, I37lX*2( 1PE11.3) , IX, 313) ,/) 

6000 FORMAT< 10X, ’***WARNING: INCONSISTENT UNDEFINED SYMMETRY’ , 

* ’ FOR PLANE’ , 12) 

7000 FORMAT! //,5X, ’ANALYTICAL VS. EXPERIMENTAL SYMMETRY: ’ ) 

8000 FORMAT( 1H1 , * NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 

* ’CORRELATION’, /IX, A9) 

END 
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SUBROUTINE INTERSECT( AVECT, EVECT, IDDOFA, IDDOFE, ICOMPA, ICOMPE, 

* NDOFA , NDOFE , I DDOF I , I COMP I , NDOF I , 

* ASYMM, ESYMM, LSYMM) 

COPY INTERSECTION SET OF VECTORS IDDOFA 8 IDDOFE TO VECTOR 
I DDOF I. NDOF I IS THE LENGTH OF THE INTERSECTION SET. 

REDUCE AVECT 8 EVECT VECTORS TO CONTAIN ONLY THE CORRESPONDING 
INTERSECTION SET OF TERMS, I.A.W. THE MASTER SET IDDOFI. 

IF THE INTERSECTION SET IS NULL (NO EQUIVALENT D.O.F ID’S), 

THE VALUE RETURNED FOR NDOF I IS 0 . 

IF SYMMETRY IS TO BE CONSIDERED, AND ANALYTICAL 8 EXPERIMENTAL 
SYMMETRY ARE NOT THE SAME, THE VALUE RETURNED FOR NDOF I IS 0 . 

REAL AVECTC *) , EVECT( *) 

INTEGER IDDOFA(*) , IDDOFE(*) , IDDOFH *) 

INTEGER ICOMPAC*) , ICOMPE( *) , ICOMPH*) 

INTEGER ASYMMC3) ,ESYMM(3) 

LOGICAL LSYMM 

IA = 1 
IE = 1 
NDOF I = 0 
IF (LSYMM) THEN 

SYMMETRY IS TO BE CONSIDERED: 

IF A 8 E ARE UNLIKE, LEAVE INTERSECTION SET AS NULL 
IF (ASYMM(l) .NE. ESYMM(l)) RETURN 
IF ( AS YMN( 2 ) .NE. ESYMM( 2) ) RETURN 
IF ( ASYMM( 3) .NE. ESYMM( 3) ) RETURN 
END IF 

DO 100 I =1, NDOF A+ NDOFE 

IF (( IA.GT. NDOFA) .OR. ( IE. GT. NDOFE) ) GOTO 150 
IF ( ( IDDOFA( IA) .EQ. IDBOFE( IE) ) .AND. 

* ( ICOMPA( I A) .EQ. ICOMPE( IE) ) ) THEN 

NDOF I = NDOF I + 1 

IDDOFI (NDOF I) = I DDOF A( I A) 

ICOMP I ( NDOF I ) = ICONPA( I A) 

AVECT ( NDOF I ) = AVECT( I A) 

EVECT(NDGFI) = EVECT( IE) 


IA ■ IA + 1 
IE = IE +1 
ELSE IF ( IDDOFA( I A) 

.LT. 

IDDOFE( IE)) 

THEN 

IA - IA + 1 
ELSE IF ( IDDOFEC IE) 

.LT. 

IDDOFA( IA)) 

THEN 

IE = IE + 1 





END IF 

100 CONTINUE 

150 RETURN 
END 
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SUBROUTINE MATCH* CORREL , NAN, MEX, NAF IT, MEF IT) 

C 

COMMON /LIMITS/ MAXMODE , MAXDOF 

DIMENSION CORREL* MAXMODE, MAXMODE) , NAF IT( MAXMODE) , MEF IT* MAXMODE) 
C 

C DETERMINE BEST EXPERIMENTAL MODE-SHAPE MATCH (HIGHEST 
C CORRELATION COEFFICIENT) FOR EACH ANALYTICAL MODE-SHAPE. 

C 

DO 200 N= 1 , NAN 
BEST = -1.0 
DO 100 M= 1 , MEX 

IF ( ABS* CORREL* N,M)) .GT. BEST) THEN 
BEST - ABS* CORREL* N.M)) 

NAF IT* N) = M 
END IF 

100 CONTINUE 

200 CONTINUE 

C 

C DETERMINE BEST ANALYTICAL MODE-SHAPE MATCH (HIGHEST 
C CORRELATION COEFFICIENT) FOR EACH EXPERIMENTAL MODE-SHAPE. 

C 

DO 400 M= 1 , MEX 
BEST = -1.0 
DO 300 N= 1 , NAN 

IF < ABS* CORREL* N,M>) . GT. BEST) THEN 
BEST = ABS* CORREL* N.M)) 

MEF IT* M) = N 
END IF 

300 CONTINUE 

400 CONTINUE 

RETURN 
END 
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G 

C 

G 


C 

C 


C 

G 


C 


C 

1000 

1500 

2000 

3000 

3500 

4000 


4500 


5000 


SUBROUTINE OPTIONS( LSYMM, PRTAN, PRTEX, RTHRESH, LPRINT) 

DETERMINE SYMMETRY, MODE-SHAPE DUMP, S THRESHOLD OPTIONS 

LOGICAL LSYMM, PRTAN, PRTEX, LPRINT 
CHARACTER* 1 YESNO 

DETERMINE IF SYMMETRY IS TO BE CONSIDERED 
TYPE 1000 

ACCEPT 2000, YESNO 

IF (YESNO .EQ. ’Y’ ) LSYMM = . TRUE, 

IF (LPRINT) PRINT 1500, YESNO 
DETERMINE IF ANALYTICAL AND/OR EXPERIMENTAL MODE-SHAPE 

VECTORS ARE TO BE PRINTED 

TYPE *, ’ ’ ?SKIP A LINE 
IF (LPRINT) PRINT *, ’ 1 
TYPE 300O, * ANALYTICAL* 

ACCEPT 2000, YESNO 

IF ( YESNO . EQ. • Y* ) PRTAN = .TRUE* 

IF (LPRINT) PRINT 3500, ’ ANALYTICAL’ , YESNO 
TYPE 3000, ’EXPERIMENTAL’ 

ACCEPT 2000, YESNO 

IF (YESNO .EO. ’Y’) PRTEX ^ .TRUE. 

IF (LPRINT) PRINT 3500, ’ EXPER I MENTAL * , YESNO 
GET RELATIVE DEVIATION THRESHOLD. DEFAULT IS 5.0 % 

RTHRESH = . 05 
TYPE 4000, RTHRESH 
ACCEPT *, RTHRESH 
TYPE 5000, RTHRESH 
IF (LPRINT) THEN 

PRINT 4500, .05, RTHRESH 
PRINT 5000, RTHRESH 
END IF 

RETURN 

FORMAT( /, * IS ANALYTICAL VS. EXPERIMENTAL SYMMETRY TO BE 

* ’CONSIDERED? (YORN): ’,©) 

FORMAT( / , * IS ANALYTICAL VS. EXPERIMENTAL SYMMETRY TO BE ’ , 

* ’CONSIDERED? (YORN): ’,A1) 

FORMAT(Al) 

FORMAT( ’ PRINT \A12,’ MODE-SHAPE VECTORS? (Y OR N) s ’,S) 
FORMAT( ’ PRINT ’,A12,’ MODE-SHAPE VECTORS? (Y OR N) s ’,A1) 
FORMAT!//’ RELATIVE DEVIATIONS GREATER THAN A THRESHOLD’ 

* ’ VALUE WILL EE PRINTED. ’ , 

* /’ THE DEFAULT THRESHOLD IS ’,F6.3, 

* /’ ENTER DESIRED THRESHOLD. ■ , 11 FOR DEFAULT: ’,$) 
FORMAT!//’ RELATIVE DEVIATIONS GREATER THAN A THRESHOLD’, 

* ’ VALUE WILL BE PRINTED. ’ , 

* /’ THE DEFAULT THRESHOLD IS ’ ,F6.3, 

* /’ ENTER DESIRED THRESHOLD. ■ , " FOR DEFAULT: ’,F6.3) 
FORMAT! ’0 RELATIVE DEVIATIONS GREATER THAN ’,2PF7.2, 

* ’ % WILL BE PRINTED. ’ ) 

END 
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SUBROUT I NE PRI NT( N 1 , FREttl , N2 , FREQ2 , NDOF , CORREL , C , S , RMSA , RMSE , 

* DIFF.DIFABS, IDDOF, ICOMP, RTHRESH) 

C 

C PRINT CORRELATION INFORMATION FOR MODES N1,N2 COMPARISON 
C 

DIMENSION DIFF(*) ,DIFABS(*) , IDDOF(*) , ICOMP(*) 

DIMENSION ID( 30000) 

C 

C DETERMINE MAXIMUM DIFFERENCE 8 DIFFS > THRESHOLD 

DIFMAX = 0.0 
IDMAX = 0 
NDIFF = 0 

IF (NDOF .GT. 0) THEN 

DIFMAX = ABSCDIFF(l)) 

IDMAX = IDDOF( 1) 

DO 100 N= 1 , NDOF 

IF ( ABS( DIFF( N) ) .GT. DIFMAX) THEN 
DIFMAX = ABS( DIFFC N) ) 

IDMAX = IDDOF(N) 

END IF 

IF ( ABS(DIFF(N) ) .GT. RTHRESH) THEN 
NDIFF = NDIFF + 1 
ID( NDIFF) = N 
END IF 

100 CONTINUE 

END IF 
C 

C PRINT SUMMARY LINE FOR THIS MATCHED PAIR 

PRINT 1000, N1,FIIEQ1,N2,FREQ2, NDOF, CORREL, C,S, 

* RMSA, RMSE, DIFMAX, IDMAX 
IF (NDOF .Eft. 0) RETURN 

C 

C PRINT ALL DIFFS > THRESHOLD , DEFAULT IS 55? 

IF (NDIFF .GT. 0) THEN 

PRINT 2000, RTHRESH, ( IDDOF ( I D(N)) , ICOMP( ID( N) ) , 

* DIFF(ID(N)), N=l, NDIFF) 

PRINT 3000, ( IDDOF( ID(N) ) , ICOMP( ID( N) ) , DIFABS( ID( N) ) , 

* N s 1 .NDIFF) 

END IF 

RETURN 

C 

1000 FORMAT! ’O’ , I6,2(F15.6, 17) , 4X, F9 . 3, 2X,5( 1PE11.3) , 19) 

2000 FORMAT( 2 IX, ’RELATIVE DEVIATIONS ( XA/RMSA-XE/RMSE) > ’, 

* 2PF7.2, ’)? : (GRID ID/DEVIATION) ’ , 

* 20(/21X,5( 19, , 1 1, , 1PE10.3))) 

3000 FORMAT! 2 IX, ’ SCALED DIFFERENCES ( (XA-XEl/S ) > THRESHOLD 

* ’ (GRID ID/DIFFERENCE) ’ 

* 20(/21X,5( 19, , II, ' , 1PE10.3))) 

END 
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c 

c 

c 

100 

c 

c 

c 

200 

C 

1000 

2O00 

3000 

4000 

5000 


SUBROUTINE PRINTPHI ( NAN, NDOFA, IDDOFA, ICOMPA, DISPAN, PRTAN, 

* HEX, NDOFE, IDDOFE, ICOMPE, DISPEX, PRTEX) 

COMMON /LIMITS/ MAXMODE , MAXDOF 
COMMON /VHEN/ CDATE 

REAL DISPAN( MAXDOF , MAXMODE) , DISPEX( MAXDOF , MAXMODE) 

INTEGER IDDOFA( MAXDOF , MAXMODE) , IDDOFEC MAXDOF , MAXMODE) 

INTEGER ICOMPA( MAXDOF , MAXMODE) , I COMP E( MAXDOF , MAXMODE) 

INTEGER NDOFA( MAXMODE) , NDOFE( MAXMODE) 

LOGICAL PRTAN, PRTEX 
CHARACTER CDATE* 9 

IF (PRTAN) THEN ! PRINT ANALYTICAL MODE-SHAPE VECTORS 
PRINT HEADER 

PRINT 500O, CDATE 
PRINT 1000 
PRINT ALL VECTORS 
DO 100 1=1, NAN 

PRINT 2000, I, ( IDDOFAC J, I) , ICOMPA( J, I) ,DISPAN( J, I) , 

* J= 1 , NDOFA( I) ) 

CONTINUE 

END IF 

IF (PRTEX) THEN ! PRINT EXPERIMENTAL MODE-SHAPE VECTORS 
PRINT HEADER 

PRINT 5000, CDATE 
PRINT 3000 
PRINT ALL VECTORS 
DO 200 1=1, MEX 

PRINT 4000, I,( IDDOFE( J, I) , ICOMPE( J, I) , DISPEX( J , I) , 

* J= 1 , NDOFE( I) ) 

CONTINUE 

END IF 
RETURN 

FORMATC 1 H0 , ’ 6 . ANALYTICAL MODE-SHAPE VECTORS * , 

* ’ ( GRID PT/D I SPL ACEMENT) ’ ) 

FORMAT( 1H0 , ’ANALYTICAL MODE’ , 14, 

* 200(/5X,5( 19, ’-Ml, ’/’ , 1PE10.3))) 

FORMAT( 1H0, ’7. EXPERIMENTAL MODE-SHAPE VECTORS 

* ’ (GRID PT/D I SPLACEMENT) ’ ) 

FORMAT( 1H0, ’EXPERIMENTAL MODE’, 14, 

* 200(/5X,5( 19, , II, '/’ , 1PE10, 3) ) ) 

FORMAT( 1H1 , * NASTRAN MODAL ANALYSIS - MODAL SURVEY STATISTICAL 

* ’CORRELATION’, /IX, A9) 

END 
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SUBROUT I NE RMSD I FF ( AVECT , EVECT , NDOF I , RMS A, RMSE , CORREL , 

* S, DIFF.DIFABS) 

CALCULATE RELATIVE (TO RMS) DIFFERENCE BETWEEN ANALYTICAL S 
EXPERIMENTAL DISPLACEMENT FOR EACH GRID POINT. 

CALCULATE ABSOLUTE DIFFERENCE SCALED TO ANALYTICAL. 

REAL AVECT( *) ,EVECT(*) ,DIFF(*) ,DIFABS(*) 

IF (NDOFI .EQ. 0) RETURN 
DO 100 1=1, NDOFI 

IF (CORREL .GE. 0.0) THEN 
DIFF( I) = O.O 

IF ( (RMSA .NE. 0.0) .AND. (RMSE .NE. 0.0) ) 

* DIFF(I) = AVECT( I ) /RHSA - EVECT( I) /RMSE 
DIFABS( I) = 0.0 
IF (S .NE. 0.0) 

* DIFABS(I) = (AVECT(I) - EVECT(I)) / S 
ELSE IF (CORREL .LT. 0.0) THEN 

C CORRECT FOR 180-DEGREE PHASE SHIFT 

DIFF( I) = 0.0 

IF ( (RMSA .NE. 0.0) .AND. (RMSE .NE. 0.0) ) 

* DIFF(I) = AVECT( I ) /RMSA + EVECT( I)/RMSE 
DIFABS(I) = 0.0 

IF (S .NE. 0.0) 

* DIFABS( I) = (AVECT(I) + EVECT(I)) / S 
END IF 

100 CONTINUE 

RETURN 
END 
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SUBROUTINE SWAP(I,J) 

INTERCHANGE THE VALUES IN TWO 4-BYTE VARIABLES 

I TEMP = I 
I = J 
J = I TEMP 
RETURN 
END 
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LAMA 


The table of modal frequencies, masses, and stiffnesses output from NASTRAN are con- 
verted by processer program LAMA from table format to matrix format to satisfy the input 
needs of STATCORR. LAMA operates interactively and prompts the user for the input file 
name and the output file name. 
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tee 

200 c 

300 

400 

500 

600 

700 

800 G 
900 

1000 G 
1100 G 
1200 

1300 G 
1400 

1500 G 
1600 
1700 
1800 
1900 
2000 

2100 10O 

2200 G 
2300 

2400 G 
2500 

2600 C 
2700 

2800 G 
2900 
3000 
3100 


PROGRAM LAMATABLE 
PARAMETER (MAXEIG = 300) 

CHARACTER LAMA* 8 

REAL FREQ( MAXEIG) , GMASS< MAXEIG) , ST IFF ( MAXEIG) ,SYMM<3) 

REAL REC2(7*MAXEIG) 

DATA LUR_P , LUN_M /I 2/ If PUNCH 8 MATRIX FILE LOG UNIT *’S 

POSITION PUNCH FILE TO 1ST CARD OF RECORD 2 OF LAMA TABLE 
CALL FINDREC2C LUN_P, LAMA) 

READ RECORD 2 OF LAMA TABLE INTO ARRAY REC2 
CALL RE ADREC2 ( REC2 , RVALUES , LUR_P , MAXE I G) 

EXTRACT FREQUENCY, MASS, 8 STIFFNESS VECTORS FROM REC2 
NEIGERS = NVALUES/7 
DO 100 1=1, RE I GENS 

FREQ(I) = REC2(7*I-2) 

GMASS(I) = REC2( 7*1- 1 ) 

STIFF(I) = REC2C 7* I ) 

CONTINUE 

DO ASCENDING SORT OF FREQ VECTOR; SLAVE SORT MASS 8 STIFF VECTORS 
CALL SORTCNEIGENS, FREQ, GMASS, STIFF) 

CHECK PUNCH FILE FOR SYMMETRY PARAMETERS 
CALL READS YM( SYMM, LUN_P) 

CLOSE PUNCH FILE 

CLOSE ( UN I T= LUN_P ) 

WRITE MATRIX FILE 

CALL WRITEMTXC NEIGENS , FREQ, GMASS , STIFF , SYMM, LUN_M, LAMA) 

STOP 

END 
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100 

200 

C 

300 

400 

C 

500 

C 

660 

10 

700 

800 

900 

c 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

100 

1700 

C 

1800 

125 

1900 

2000 

2100 

C 

2200 

150 

2300 

2400 

2500 

2600 

200 

2700 

C 

2800 

225 

2900 

3000 

3100 

C 

3200 

250 

3300 

3400 

C 

3500 

1000 

3600 

2000 

3700 

3000 

3800 

4000 

3900 



SUBROUTINE F I NDREC2 C LUN—P , LAMA) 
CHARACTER PF I LE*80 , LAMA*8 , L I NE*80 


OPEM -mE C ?«0^ E ’ ENTER I SPOT POUCH PILE HAME! * 

ipST5«?Kj« HI mME-PFILE, TYPE- • OLD’ , READOHLY. ERR= 10) 

P0S1 ^ Uw.^E^KoF LAMA TABLE: ■ 

ACCEPT 3000, LAMA 

DO 100 1=1,1000000 . ,„ F 

“*?! tSif ) «D. < LIHEC 0 ! I«> ! EO.LAMA. . 

* * GOTO 150* n LAMA TABLE FOUND 

CONTINUE ~ 

TUU NOT FOUND' 

CLOSE! UN IT=LUH_P) 

POSITION^C^FILE TO RECORD 2 OF LAMA TABLE 

. 00 

CONTINUE 

: 0T .SS D k S C™ T 2 °™LAMA TABLE HOT FOOHD’ 
CLOSE! UNI T=LUN_P> 

POSITION FILE TO REREAD 1ST CARD OF RECORD 2 
BACKSPACE LUN-P 

BF.TIIRN 


FORMAT! IX, A30, IX,*) 
FORMAT! A80) 

FORMAT! A8) 

FORMAT! A24, 116) 

END 
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100 

200 

300 

400 
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600 
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2800 
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3300 

3400 


C 

C 


100 

200 

C 

225 

C 

250 

C 

1000 

2000 

3000 

4000 


SUBROUTINE READREC2( IREC2, RVALUES, LUN_P, MAXEIG) 

INTEGER IREC2!7*MAXEIG) , ITEMPC4) 

CHARACTER* 6 CTEMP(4) 

READ J,®X CARD 0F LAMA TABLE RECORD 2 

,HEC2 ‘ 

RE ' U, DfMo"!”?.Ser r,L F0 ™ D 

Do''?e0™M , 4 OO0 ’ E " D ' 225> <CTE,IP< J > • itempc J) , j, 

iWssr&A ; E ^ ,Ec ' , 0010250 

I REC2(N VALUES) = ITEMP(J) 

IF (CTEMP( J)(6:6) .Ett. *-*) 

I REC2 C NVALUES ) = — IREC2! NVAI TJF<31 
CONTINUE 1 ^ 8 * EQ * 7 * MAXEIG) GOTO 225 
CONTINUE 

END tWSSS^hSlSS 4 table not found 

CHFCK’ Tn P c?r E ,r TYPE "CONTINUE" OR "STOP"* 

return -SnMt^nT BY ’ 

FORMAT! 40X, 2116) 

FORMAT! 8X, 4! A6 , 110)) 

^FORMATUX, !*** WARNING: END OF RECORD 2 NOT REACHED.*, 

FORMAT! IX,’ ’***VARNING: # WMCORD^^ENTRIE^ ‘ 14 ’ I4> 
END N ° T DIVISIBLE BY 7*) " IES ,14, 


1,4) 
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100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 
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1300 

1400 

1500 

1600 

1700 

1800 

1900 

2O0O 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

34O0 

3500 

3600 

3700 

3800 


SUBROUTINE READS YM( SYMM, LUN_P> 

C 

CHARACTER* 80 RECORD 
REAL SYMM(3) 

LOGICAL LSYM1 , LSYM2 , LSYM3 
C 

REWIND LUN_P ! START SEARCH FROM TOP OF FILE 
C SET SYMMETRY DEFAULTS 
SYMM( 1 ) = 0.0 
S YMM( 2 ) = 0.0 
SYMM( 3) = 0.0 

C SEARCH PUNCH FILE FOR SYMMETRY DEFINITION RECORDS 

DO 10© 1=1,1000000 

READ( LUN-P , 1000, END= 150) RECORD 
IF ( RECORD( 2 1 : 28) . EQ. ’ SYH1PLAN’ ) THEN 
LSYM1 = .TRUE. 

DECODE( 2 , 2000 , RECORD( 47 s 48) ) SYMM( 1 ) 

TYPE *, ’PLANE 1 SYMMETRY:’, SYMM( 1) 

ELSE IF ( RECORD( 2 1 : 28) .EQ. ’ S YM2PLAN ’ ) THEN 
LSYII2 = .TRUE. 

DECODE( 2 , 2000 , RECORD ( 47: 48) ) SYMM( 2) 

TYPE *, ’PLANE 2 SYMMETRY:’, SYMM(2) 

ELSE IF ( RECORD* 2 1 : 28) . EQ. ’ SYM3PLAN’ ) THE:! 
LSYM3 = .TRUE. 

DECODE( 2 , 2O00 , RECORD( 47 : 48) ) SYMM( 3) 

TYPE *, ’PLANE 3 SYMMETRY:’, SYMM(3) 

END IF 


100 

CONTINUE 



c 

REPORT NON-DEFINED SYMMETRY 


150 

IF (.NOT. LSYM1) 

TYPE 3000, 

’ 1’ 


IF (.NOT. LSYM2) 

TYPE 3000, 

’2’ 

C 

1000 

2000 

IF (.NOT. LSYM3) 
RETURN 

FORMAT( ABO) 

FORMAT ( F2.0) 

TYPE 3000, 

’3’ 

3000 

FORMAT( IX, ’ PLANE ’ , A1 
END 

, ’ SYMMETRY 

NOT SPECIFIED’) 
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100 


SUBROUTINE SORT( N , FREQ, GMASS f STIFF) 

200 

C 




300 


REAL FRECKN) ,GMASS(N) # STIFF! N) 


400 

C 




500 

C 

PERFORM ASCENDING SORT OF FREQ 

VECTOR 

600 

c 

SLAVE ! 

SORT GMASS AND STIFF VECTORS 

700 

c 




800 


IF (N 

.EQ. 1) RETURN 


900 

c 




1000 


DO 100 

I=1,N-1 


1100 


DO 100 

K= I + 1 1 N 


1200 


IF 

(FREQCIO .LT. FREQ(I)) 

THEN 

1300 



TEMP = FREQ< I ) 


1400 



FREQ(I) = FREQ(K) 


1500 



FREQ( IO = TEMP 


1600 



TEMP = GMASS(I) 


1700 



GMASS(I) = GMASS ( K) 


1800 



GMASS (K) = TEMP 


1900 



TEMP = STIFF(I) 


2000 



STIFF(I) = STIFF(K) 


2100 



STIFF(K) = TEMP 


2200 



END IF 


2300 

100 

CONTINUE 


2400 


RETURN 



2500 


END 
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lee 

200 C 

300 

400 

500 C 
600 C 
700 
800 
900 
1000 

1100 G 
1200 
1300 
1400 

1500 C 
1600 
1700 
1800 
1900 
2000 
2100 
2200 
2300 

2400 C 
2500 1000 

2600 2000 
2700 3000 

2800 4000 

2900 


SUBROUTINE WRITEMTX! N , FREO, GMASS , STIFF , SYMM, LUN-M, LAMA) 

CHARACTER MFILE*40,LAMA*8 

REAL FREQ(N) ,GMASS(N) , STIFF! N) , SYMM! 3) 

OPEN MATRIX FILE , 

TYPE 1000, ’ENTER OUTPUT MATRIX FILENAME:’ 

ACCEPT 2000, MFILE 

OP EN< UN I T= LUN-M , N AME= MF I LE , TYPE= ’ NEW’ , 

* CARRI AGECONTROL= ’LIST’ ) 

WRITE MATRIX FILE HEADER 
IROWS = N 
JCOLS = 6 

WRITEC LUN_M, 3000) LAMA, IROWS, JCOLS 
WRITE VECTORS AS THREE COLUMNS OF MATRIX 
WRITE! LUN-M, 4000) FREQ !! COLUMN 1 

WRITE! LUN-M, 4000) GMASS II COLUMN 2 

WRITE! LUN-M, 4000) STIFF II COLUMN 3 

WRITE! LUN-M, 4000) ( SYMM! 1) , 1= 1 , N) l I 
WRITE! LUN-M, 4000) ! SYMM! 2) , 1= 1 ,N) I! 

WRITE! LUN-M, 4000) (SYMM! 3) , 1= 1 ,N) l! 

CLOSE! UN IT= LUN-M) 

RETURN 


jOLUMN 4 
COLUMN 5 
COLUMN 6 


FORMAT! 1X.A30, IX, «) 
FORMAT! A40) 

FORMAT! A8, 2 18) 
FORMAT! 1 PE 16, 8) 

END 
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UNPACKDMI 


The data blocks of mode shapes and component degrees of freedom as output from 
ASTRAN are in a special format called DMI. This format is not acceptable to STATCORR, 
so these data must be reformatted. The pre-processor program UNPACKDMI is used to take 
modes and freedoms out of DMI format. The output of the modes from this program will be 
in final form for STATCORR. The output of freedoms will need further processing by another 
pre-processor. UNPACKDMI operates interactively and prompts the user for responses. 
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4000 
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42O0 

4300 

4400 
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4660 

4700 

4800 

4900 

5O00 

5100 

5200 

5300 

5400 

5500 

5600 

5700 

5800 

5900 

6000 

6100 


C 

C 


C 

C 


C 

399 


4O0 

C 

450 

C 

499 


5m 

550 

C 

1000 

2000 

3000 

4000 

5000 

6000 


') ALL = .TRUE. 


PROGRAM UNPACKDMI 
PARAMETER ( MAXMTXS IZ = 100000) 

CHARACTER MATNAM* 8, FLAG* I , C0NTID*3 
REAL ARRAY( MAXMTXS I Z) 

LOGICAL ALL, EOF 
DATA LUN / 10 / 

DO 500 NF I LE= 1 , 1000000 
GET NEXT FILENAME 

CALL GETFILC LUN, *499) 

ALL = .FALSE. 

EOF = .FALSE. 

DO 400 NMATRX=1, 1000000 
GET NEXT MATRIX NAME 

IF (.NOT. ALL) THEN 
TYPE 1000 

ACCEPT 2000, MATNAM 
IF (MATNAM . EQ. * *ALL 
END IF 

SEARCH FILE FOR MATRIX 8 GET DIMENSIONS 

I F L ('E0F? m WTO 1 399 ’ ’ NC0LS * C0NTID • LOU, ALL, EOF) 

IF (MROWS*NCOLS .CT. MAXTTTXSIZ) THEN 

TYPE *, ’MATRIX IS TOO LARGE’, MATNAM, MROWS, NCOLS 

*< INCREASE PARAMETER MAXMTXS IZ IN UNPACKDMI 1 
GOTO 399 
END IF 

UNPACK MATRIX INTO ARRAY 

CALL UNPACK! ARRAY, MROWS, NCOLS, CONTI D LUN ALL FOF) 
REDUCE SQUARE DIAGONAL MATRIX TO VECTOR ? ’ 

IF (MROWS .EQ. NCOLS) THEN 
TYPE 6000, MATNAM 
ACCEPT 4000, FLAG 
IF (FLAG .EQ. ’Y’) 

* „„„ CALL REDUCE! MATNAM, ARRAY, MROWS , NCOLS) 

END IF 

WRITE UNPACKED MATRIX TO NEW FILE 

CALL WRTMTX( MATNAM, ARRAY, MROWS , NCOLS) 

IS ANOTHER MATRIX IN THIS FILE DESIRED ? 

IF (ALL .AND. EOF) COTO 450 
IF (.NOT. ALL) THEN 
TYPE 3000 
ACCEPT 4000, FLAG 
IF (FLAG .NE. ’ Y’ ) GOTO 450 
END IF 

CONTINUE 

NO MORE MATRICES DESIRED FROM THIS FILE 
CLOSE (UNIT=LUN) 

IS ANOTHER FILE DESIRED ? 

TYPE 5000 
ACCEPT 4000 , FLAG 
IF (FLAG .NE. ’Y’ ) GOTO 550 
CONTINUE 

STOP ’NO MORE FILES REQUESTED.’ 

FORMAT! AH)"^ 7 *** 1 rlATRIX NAME ’ ( "* AL L’ FOR ALL MATRICES): ’ ,S) 

F0RMAT( A 1 ?° Y ° U WMT AH0THER MATRIX IN THIS FILE ? (Y OR N) : ’,8) 

F0RMAT( ’ DO YOU WANT ANOTHER FILE ? (Y OR N) : ’8) 

F0RMAT( ’ REDUCE DIAG MATRIX ’ ,A8, ’ TO A VECTOR ? (Y OR N) : ’8) 

END 
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C 


C 

C 


100 

150 

C 

C 

199 


C 

1000 

2000 


SUBROUTINE FNDMTX( MATNAM, MROWS , NCOLS , CONTID , LUN * ALL ♦ EOF) 
CHARACTER MATNAM*8, RECORD*80, FLD1T03*24, CONTID*3 
LOGICAL ALL, EOF 


IF (.NOT. ALL) THEN 
REWIND LUN 
EOF = .FALSE. 

END IF 

SEARCH FOR DESIRED MATRIX 8 GET DIMENSIONS 
FLD1T03 = * DMI f //MATNAM// * 0* 

DO 100 1=1,1000000 

READ( LUN , 1000 , END= 199 ) RECORD 
IF (ALL) FLD1T03( 9 t 16) = RECORD! 9: 16) 

IF (RECORD! 1:24) .EQ. FLD1T03) THEN 
IF (ALL) MATNAM = RECORD(9:16) 

READ( RECORD, 2000) MROWS, NCOLS, CONTID 
CONTID! 1:1) = ’*’ 

GOTO 150 
END IF 
CONTINUE 
RETURN 


END OF FILE; NEXT MATRIX NOT FOUND 
EOF = . TRUE . 

IF (ALL) TYPE *, ’END OF FILE REACHED* 

IF (.NOT. ALL) TYPE *, MATNAM, ’MATRIX HEADER 
RETURN 


CARD NOT FOUND’ 


FORMAT! A80) 

FORMAT! 56X, 218, A3) 
END 
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100 


SUBROUTINE GETFILI LUN, *) 


200 


CHARACTER FILNAM*40 ■ ... 

i 

300 

C 



400 

C 

GET NEXT FILENAME 


500 


TYPE 1000 


600 


ACCEPT 2000, FILNAM 


700 

C 

OPEN FILE 


800 

900 


OPEN! UNIT=LUN,NAME=FILNAM,TYPE= ’OLD’ .READONLY, ERR= 199) 
RETURN 


1000 

C 



1 100 

C 

NEW DMI FILE NOT FOUND 


1200 

199 

TYPE *, FILNAM, ’FILE NOT FOUND’ 

- 

1300 


RETURN 1 

- 

1400 

C 

‘ 

- 

1500 

1000 

FORMATS * ENTER NAME OF NEXT DMI FILE: ’3) 


1600 

2000 

FORMATC A40) 


1700 


END 

= 
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C 

C 


C 

c 


100 


SUBROUTINE REDUCE! MATNAM, MATRIX, MR0W3 , NCOLS) 

CHARACTERS 8 MATNAM 

REAL MATRIX* MROWS, NCOLS) 


CHECK FOR SQUARE MATRIX 

IF (NROWS .NE. NCOLS) THEN 

TYPE *, MATNAM, ’NOT A SQUARE MATRIX. 

RETURN 

END IF 


NO REDUCTION DONE.’ 


REDUCE TO (IX NCOL) 
DO 1O0 N=l, NCOLS 
MATRIX! N, 1) = 
CONTINUE 
MROWS = 1 
RETURN 

END 


ROW VECTOR USING DIAGONAL VALUES 
MATRIX! N,N) 
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lee 

200 

300 

400 

300 

600 

700 

800 

900 

1000 

1100 

1200 

1300 G 
1400 C 
1300 
1600 
1700 

1800 100 
1900 C 
2000 G 
2100 
2200 
2300 

2400 G 
2300 
2600 
2700 

2800 G 

2900 

3000 

3100 

3200 

3300 G 

3400 C 

3500 

3600 

3700 

3800 

3900 

4000 

4100 

4200 

4300 200 

4400 C 
4500 

4600 G 

4700 20 1 

4800 

4900 

500O 

5100 

5200 

5300 

5400 

3500 

5600 300 

5700 

5800 G 
5900 
6000 
6100 


SSSU*"" 1 ' imovs.BcoLS.coirrrD, um, all, eon 

CHARACTER* 1 6 FLDID(2:5), CARDID*8, C0NTID*3 
EQUIVALENCE ( RECORD! Is 1), CARDID(lsl) ) 

EQUIVALENCE ( RECORD! 9: 9), FLDID! 2) (1:1) ) 

CHARACTER* 6 FMT! 6) , REAL, INTGR, BLANK 
REAL MATRIX! MROWS, NCOLS) 

LOGICAL ALL, EOF 


DIMENSION VALUE! 4) , IVALUE(4> 

EQUIVALENCE (VALUE! 1), IVALUE(l)) 

DATA FMT! 1)/’ (8X, ’/ , FMT(6)/’8X> ’/ 

DATA REAL/ ’ E 1 6 . 8 , ’ / , INTGR/’ I 16, ’/ , BLANK/’ A16 ,’ / 


ZERO THE MATRIX 

DO 100 1=1, MROWS 
DO 100 J=l, NCOLS 

MATRIX! I, J) =0,0 
CONTINUE 


PROCESS THE DMI CARDS FOR THIS MATRIX 
DO 400 K= 1,1000000 

READ! LUN, 1000, END= 425) CARDID, FLDID 
IF (CARDID ,EQ. 'DMI ’) THEN 

TH!S IS HEADER CARD OF NEXT MATRIX; THIS MATRIX DONE. 
IF (ALL) BACKSPACE- LUN 
GOTO 450 

ELSE IF (CARDID , EQ, ’DMI* ’) THEN 
THIS IS A COLUMN IDENTIFIER CARD 
READ! RECORD, 2000) JCOL, IRON, VALUE! 1 ) 

MATRIX! IRON, JCOL) = VALUE! 1) 

IRON = IRON + 1 

ELSE IF (CARDID! 1:3) .EQ, CONTID) THEN 
THIS IS A RON/VALUE CARD 
DETERMINE TYPE OF CONTENTS IN FIELDS 2:5 
DO 200 L=2, 5 

IF ( FLDID! L)( 13: 13) . EQ. ’E’) THEN 
FMT! L) = REAL 

ELSE IF ( FLDID! L) ( 16 : 16) .EQ. ’ ’) THEN 
FMT! L) = BLANK 

ELSE 


FMT! L) = INTGR 
END IF 
CONTINUE 

READ CARD IN APPROPRIATE FORMAT 

READ! RECORD , FMT , ERR= 20 1 ) VALUE 
PROCESS FIELDS 2:5 
DO 300 L=2, 5 

IF (FMT(L) .EQ. REAL) THEN 

MATRIX! IRON, JCOL) = VALUE! L- 1 ) 
IRON = IRON + 1 

ELSE IF (FMT(L) .EQ. INTGR) THEN 
IRON = IVALUE(L-l) 

ELSE IF (FMT(L) .EQ. BLANK) THEN 
CONTINUE 
END IF 
CONTINUE 

ELSE 


THIS IS A NON- DMI CARD; THIS MATRIX DONE 
READ! RECORD, 1000 ) CARDID, FLDID 
GOTO 450 
END IF 
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6200 

400 

6300 

C 

6400 

425 

6500 

C 

6600 

450 

6700 

C 

6800 

1000 

6900 

7000 

2000 


CONTINUE 
END OF FILE 

EOF = .TRUE. 

END OF MATRIX 
RETURN 

FORMAT! A8.4A16) 

FORMAT! 24X.2I 16.E16.8) 
END 


100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 


C 


C 

C 


C 

1000 

2000 


SUBROUTINE WRTMTXC MATNAM, MATRIX, MROVS, NCOLS) 
CHARACTER MATNAM*8, NAME* 12 
REAL MATR I XC MROWS , NCOLS ) 

NAME = MATNAM// ’ . MTX' 

NA ME=NAME, TYPE= ’ NEW , CARR I ACEC0NTR0L 3 
WRITE HEADER RECORD TO ASCII FILE 

WRITE! 50, 1000) MATNAM, MROWS, NCOLS 
WRIT L M ^I RIX TO ASCII FILE ( 1 VALUE PER RECORD ) 
WRITE! 50, 2000) MATRIX 
CLOSE !UNIT=50) 

TYPE *, ’MATRIX FILE WRITTEN WITH MROWS, NCOLS:’ 

TYPE *, NAME, MROWS, NCOLS 

RETURN 


FORMAT! A8.2I8) 
FORMAT! 1 PE 16. 8) 
END 


’LIST’) 
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GRDPTLST 


The freedoms output from processor UNPACKDMI are input to processor GRDPTLST to 
put them into final format for STATCORR. GRDPTLST operates interactively and prompts 
the user for the input file name and the output file name. 
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PROGRAM GRDPTLST 

C THIS PROGRAM WILL READ THE TESET (DOF ID VS. GRID POINT ID) 

C MATRIX FILE OBTAINED FROM THE NASTRAN PUNCH FILE AND GENERATE 

C A GRID POINT ID LIST FILE, USING THE NON- ZERO ENTRIES OF THE 

C TESET MATRIX. 

C 

CHARACTER*35 MTXFIL, GPFIL, ITTXNAME*8 


C 

C 

10 


C 


C 


100 


C 

1000 

2000 

3000 

4000 

5000 

6000 


GET TESET MATRIX FILENAME 

TYPE 1000, * ENTER GRID POINT ID MATRIX FILENAME:’ 

ACCEPT 2000, MTXFIL 

OPEN( UNIT= 1 , NAME=MTXFIL,TYPE= ’OLD’ , READONLY, ERR= 10) 

READ( 1,3000) MTXNAME, MROWS, NCOLS 

IF (NCOLS .NE. 1) STOP ’ERROR: SHOULD BE ONLY ONE COLUMN’ 
GET GRID POINT LIST FILENAME 

TYPE 1000, ’ENTER GRID POINT ID OUTPUT LIST FILENAME:’ 
ACCEPT 2000, GPFIL , 

0PEN( UNIT=2, NAME=GPFIL, TYPE= ’ NEW’ , CARR I AGECONTROL= ’ LIST’ ) 
GENERATE LIST FILE FROM NON- ZERO MATRIX ENTRIES 
NIDS = 0 

DO 10O MDOF= 1 , MROWS 

READC 1 , 5000) GRDPT 
IF (GRDPT .NE. 0.0) THEN 
I GRDPT = INT( GRDPT) 

COMP = GRDPT - FLOAT( I GRDPT) 

ICOMP * NINTC 10.0*COMP) 

IF (ICOMP .GT. 6) THEN 
I GRDPT * I GRDPT + 1 
ICOMP = O 
END IF 

WRITE( 2 , 6000) I GRDPT, ICOMP, MDOF 
NIDS = NIDS + 1 
END IF 
CONTINUE 
CLOSE( UNIT r 1) 

CLOSE( UN I T= 2 ) 

TYPE *, NIDS, ’ENTRIES WRITTEN TO GRID POINT LIST’ 

STOP 


FORMATC IX, A44, IX, 6) 

FORMATC A35) 

FORMATC AS, 218) 

FORMAT< IX, A8, ’ IS A ’,18,’ ROW BY ’,18,’ COLUMN MATRIX ) 
FORMATC E16. 8) 

FORMATC 18, ’-’,11, 112) 

END 
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FUTURE REQUIREMENTS 


1 . If the practice of comparing analysis with test, as outlined here, for certifying a design 
should become commonplace, one could imagine that the demand for additional information 
would increase. Instead of depending on a comparative analysis to indicate that differences do ex- 
ist, and roughly where, engineers increasingly might require comparative analysis to predict more 
precisely where differences arise and their causes. An extension of the existing method in this di- 
rection would be to apply a correlation function using a correlation displacement variable. Con- 
sider the expression 


0 (P,r,s) 



T(x p ,y p ,z) x A[(x p 


r, y p - s)z] dV. 


0 measures the best correlation of a test function T at a particular point P(x p ,y p ,z p ) with an analy- 
tical function A in the neighborhood of P:[V p ] as the analytical function is evaluated at displace- 
ment distances (r and s) away from P in x and y. The motivation for considering such a quantity is 
to explore a point at which earlier results of a correlation coefficient showed consistent differences 
over a number of modes at one location of the structure. Since tests are limited as to the intensity 
of data obtainable in a small region, the test coordinate can be held constant in one or two dimen- 
sions while the analysis mesh is intensified in the region to provide enough definition to step 
through values of r and s to find what location in the vicinity of P shows greater correlation of A 
with T at P. Such a spatial correlation function might prove helpful in exploring a design defect 
or an analytical modeling error of a structure. 

2. Another extension is to equip the correlation to make comparisons on the basis of mode 
shapes defined in the complex domain. Correlation with a damping displacement variable could 
prove useful in refining the damping away from average modal considerations to spatial distribu- 
tions of damping. 

3. Currently, the analysis presumes that any particular point on the test article, instrumen- 
tation can detect only one component direction of behavior. If multi-axis behavior at a test point 
can be provided, the component at an analytical point can be implemented by adding one more 
digit after the decimal point of the real number representing a grid point. The processor program 
can discriminate to the nearest integer by multiplying by 10, and then rounding off. Thus, the 
real value in TESET and DOFLIST can reflect correspondence in keeping with increased experi- 
mental detail. 

4. Processing of analytical data in preparation for STATCORR is handled in discrete op- 
erations. If a number of different analytical runs are being processed at the same time for later 
assembly as STATCORR input, their separate results might get uncoordinated, so that a symmetric 
mode shape, adjusted for sign at certain image points might be identified with a wrong frequency 
or modal stiffness. With further effort, the several operations could be combined with automated 
bookkeeping to eliminate the risk of getting results scrambled. 

?-■ Human error can be further reduced if the entire analysis were made interactive with 
computer prompts. Then when STATDMAP is accessed from USERLIB, the whole sequence of 
orderly information train could be set in motion to keep all relations coordinated. 


84 



6. STATCORR can be readily extended to comparisons of static loadings. First of all, the 
mean, x', has entirely different meaning in statics than in modal analysis, x' for statics cannot be 
prescribed as the equilibrium position. The non-zero mean must be reintroduced into the sto- 
chastic quantities. Scaling was meaningful for eigenvalues because the method of normalizing 
could be arbitrary. The relative deviation of equation (7) should use actual unsealed response data 
for identifying magnitudes of difference. 

7. The implementation of the option for comparing symmetrically based analysis against 
samples over the entire test article was restricted to the use of basic coordinates only. A logical 
extension is to allow any number of rectangular coordinates for defining direct and image points. 
Enough extra work would be involved to implement more than one coordinate system, that it 
may be tolerable to limit the First extension of rectangular coordinates whose planes are parallel 
to the planes of symmetry. In a second extension, all types of coordinates can be logically ad- 
mitted. 


USER MANUAL SUMMARY 


This document is intended to serve as a guide for comparing analytical modes with test 
when using NASTRAN, STATCORR and its associated processors, in conjunction with the test 
data. 

First, select one of three bases of comparison: (1) whole vs. whole in which both the analyti- 
cal modeling and the test instrumentation involves the whole structure; (2) symmetric vs. sym- 
metric in which the analytical modeling and the test instrumentation involves only a basic segment 
of reflective symmetry; and (3) symmetric vs. whole in which the analysis is based upon modeling 
of only a basic segment of reflective symmetry and the test is based upon instrumenting the 
whole structure. 
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USER MANUAL SUMMARY 
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